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INTRODUCTION 


The United States currently utilizes a rather small family of launch 
vehicles (boosters) to support a \ '.ried spectrum of satellite and spacecraft 
programs tl]. These launch vehicles have been carefully designed to 
accommodate a wide range of payload configurations. In general, the payload 
interfaces with the launch vehicle at a limited subset of candidate structural 
”hard** points at the payload launch vehicle separation plane. The latest 
example in the series is the STS (Space Transportation System). 

It is important that any candidate payload be designed to withstand the 
load environment transmitted to the payload from within the shielded payload 
compartment* Such environments coimaonly originate from a static (steady 
state) vehicle acceleration, a transient or dynamic event such as rocket motor 
ignition, or an acoustical environment. Very often, it is the transient 
dynamic response behavior of the payload that constitutes payload design load 
profiles; hence, it is important that proper attention be given to the payload 
transient response characteristics as influencing major design decisions. As 
an example, let us consider the landing of the orbiter (= a delta-wing- 
airplane-like module) carrying a certain payload. Obviously, when the orbiter 
touches the ground it will experience reaction forces. These forces will be 
transmitted to the payload through the interface (i.e. through the connection 
points between the orbiter and its payload). The payload will undergo 

elastic displacements. The question is, will the payload be able to withstand 
those displacements without being damaged? The answer to this question 
requires a dynamic analysis of the booster/payload system as we shall see in 
subsequent sections. This dynamic analysis will yield the elastic 
displacements in the payload due to the known reaction forces on the orbiter. 
These displacements can then be used to calculate the internal forces in the 
different members of the payload which in turn leads to the calculation of 
stresses and strains in those members. Finally, these stresses and strains 
enable the payload designer to determine whether or not the members will be 
damaged during the landing event. 
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LJ’resent analytical techniques by which such design loads are predicted are 
very costly and time consuming^ A typical cycle (as the above mentioned 

process is called) generally requires: 

1* Generation of a payload model; 

2. Calculation of the modal characteristics of the payload restrained at 
the interface; 

3* Formation of a transformation to couple the payload to the booster 
interface; 

4. Coupling of the payload to the booster and calculation of the system 
modal characteristics; 

5. Calculation of the time response of Lhe system to the specified 
forces; 

6. Use of the time response results to calculate loads. 

The calendar turnarouuu time of a given cycle usually is lengthened when 
the payload design organization, the booster organization and the payload 
integration organization are different companies. The reason for this is that 
a fair amount of coordination is necessary to make the transfer of information 
between those three organizations optimal. Unfortunately, this coordination 
is very difficult to establish resulting in considerable time delays. 

Moreover, these costs and delays repeat themselves for every load cycle (i.e., 
every time a change is made in the booster or payload). A typical example is 
the development of the Viking Orbiter System [2]. Upward of nine 
organizations were responsible for hardware or integration functions which 
directly affected the evaluation of dynamic transient loads. The number of 
interfaces between those organizations resulted in difficulties in arranging 
for the necessary analyses at each organization, in obtaining the necessary 
data, in establishing priorities, in establishing output requirements, and in 
correctly transferring data between organizations. The time duration for one 
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load cycle ranged from three to twelve months which depended on the nund)er of 
events, forcing functions per event, and complexities of the analysis. Of 
course, if the booster already has its final design, many of these problems 
can be avoided. Theoretically, only one transfer of booster data to the 
payload organization would be necessary. 

It is clear that the need exists for a so-called ** short-cut** methodology 
in this area, A **short-cut** method should meet three essential requirements. 

1, It should take advantage of the fact that the booster stays the same from 
one design cycle to the next or from one flight to the next. The payload 
integrator should be able to reuse several previously calculated booster 
quantities (e.g, mass and stiffness matrices, modes etc.) 

2, The ** short-cut** method should avoid as much as possible the transfer of 
information between different organizations involved in the load 
analysis. Ideally, the payload organization should be able to estimate 
design loads to support their design activities without having to rely on 
other organizations, A one-time transfer of booster information should 
suffice if the payload integration organization is the same as the payload 
organization [ 2 J , 

3, A **short-cut** method should be cost-effective* For example, it seems 
reasonable that no complete cycle (i.e, item 1-6 on previous page) is 
necessary if only small changes are made in the payload. A similar 
situation exists in the assessment of STS payload design loads. Although 
many of these payloads will be designed for multiple flights with moderate 
changes, state of the art dynamic loads prediction technology does not 
provide a way to avoid complete reanalysis of the booster/payload system 
[ 1 ]. 

The object of this contract then is, to develop and implement such a 
**short-cut** methodology. The present Assessment Report covers Study Task I of 
the contract. Chapter I presents the standard techniques used to analyze a 
payload/booster system. They are **full-scale** methods in the sense that they 
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all require the solution of the coupled equations of motion of the 
booscer/payload system. Chapter II identifies several "short-cut” 
methodologies. These already existing techniques do not require the solution 
of the coupled system equations. The potentials and shortcomings of each of 
these methods is discussed. Chapter III covers the "favored" methods 
accompanied by recommendations for further development, refinement and 
demonstrations. We also included the outline of a new approach. 
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CHAPTER I; STANDARD TECHNIQUES 


1 . INTRODUCTION 


The objective of this first chapter is to identify and assess the most 
prominent standard techniques currently available. This will allow the 
introduction of the necessary background information in terms of a unified 
nomenclature. It will also provide us with the state-of-the-art full-scale 
methodology. Such a method is necessary for comparison purposes. Also, some 
of the features of these methods may be incorported into some of the short-cut 
methods. This chapter will give us the opportunity to more clearly identify 
the requirements of an acceptable short-cut methodology. Tne first section 
deals with the equations of motion in the discrete time domain. 

2. THE EQUATIONS OF MOTION IN THE DISCRETE TIME DOMAIN [3,4] 


The objective of this section is the derivation of the ec^pations of motion 
of the booster/payload system. Figure 1 shows the free body diagrams of the 
booster B and the payload P. The booster and the payload are connected to 
each other through the interface. Physically, the interface is the collection 
of structural "hard** points which the booster and the payload have in common. 
Mathematically, this means that the generalized displacement vector 




on the 


booster side of the interface must be equal to its equivalent | | on the payload 


side. Hence 
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for all times t 


( 1 ) 
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Figure 1 Free-Body Diagrams of Booster B and P( 'load P 
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similarly, the generalized reaction vectors|R® ^and interface satisfy, 

|r® I * _ all times t (2) 


From the free body diagrams in Figure 1 we can easily write the equations 
of motion for the booster B and the payload P as, n 


' % 




where represents the generalized displacement ^rector of B. This vector 

can be partitioned according to non-interface displacements and 

interface displacements i f 9 

' ^ ' i Tt \ 


Similarly for P, 




and the stiffness matrix 


The mass matrix stiffness 

also be partitioned in the same manner^ 


of the booster B can 


Similarly for P, 
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Finally, the vector represents the externally applied forces on the 
boooter B and can be written as 


N 


‘N 


( 8 ) 


Similarly for P, 


{^p} = 


/ pP \ 


(9) 


Using Eqs. (4), (5), (8) and (9) we can rewrite Eq. (3) as follows, 

^B 


’ ^ 




“p. 


"B 


"P 


"B 



"P 

+ 

rw 


B 

B 


P 

P 


N 


N 

' P 


( 10 ) 


Both equations (3) and (10) represent the equations of motion of the 
undamped booster B and the undamped payload P. In order to derive the 
equations of motion for the coupled system (i*e. booster/payload system), we 
need to eliminate the a priori unknown reactions | and-|k^ • We shall 
now establish a convenient and physically meaningful way to accomplish this 
elimination. To this end, let us solve the third partition of Eq. (10) for 
the non-“interface displacement vector of the payload P. 

It is now noted that the non~interface displacement vector | | consists of 

two parts. To understand the physical meaning of these two terms let us 
assume that the interface displacements are zero i.e. | *{ ^ } * that 
case it follows from Eq. (11) that the second term on the right-hand side can 
be interpreted as the non-interface displacement of the payload with respect 
to the interface. Let us denote this term by | | • It is then clear that the 

first term on the right-hand side of Eq. (11) represents tht non-interface 
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displacement of the payload due to the displacement of the interface 

only . Therefore, Eq. Ol) can be written as 


{ 4 } “ H 
[®P] ■ - [<N ] ' [>^ iJ 


Equations (12-13) are now used to establish the following transformation, 


-f 


where we used Sqs. (4), (5), (1), (12) and (13). This transformation will 
eliminate the redundant set of displacements/ >in Eq. (10) and in the 

\^; (P'l 

process it will also eliminate the unknown reactions <R^Vand2R First, Ic 
us introduce a more convenient notation, 

r I 0 |o -1 rx o| 0 ] 


0 I } 0 

-h = A 


Tp I Ip 


Olio 
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With this notation we now substitute the transformation Eq. (14) into £q* 


(10) and premultiply by A (T » transpose). This yields the following result^ 

B R 

N 

H 






-] 

‘pVp i 'J'V'p J 


/ 


p 


“i, * j 




WpJ 


N 


^ ^ s'' I-'; 

] I I* N 


V 


(17) 


mat 


At tliis point a few remarks are in order. First, we shall show that the 
rix product J ^ zero. From Eqs. (7) and (16) wc have 


W ’ H ['.1 ■ 

0 0 
rn 

r p Pi 
• Sin ‘Sii 


I 

= 



P P 

[•"in •'iiJ 


0 



.T „P 


4 


( 18 ) 


Substituting the cxnr«jOL;ion (1 3) for[ S ]ln the lower lialf oi lU;. (18) violets. 


because 


M [CJ"Kn] = -[<i] Kn] 

T \ 

l"!..] ■ Kl] (M 


[ 4 ] = [ o ] 


(19) 


Secondly, the triple matrix product 


w w w 


0 0 


“5 r.,P T r 


T 

s' I 


m Sii 


p p 

k: 


P L‘ IN ‘'II J 


0 


0 I 




( 20 ) 


will be zero for a stat::ally determinate interface. The interface is 
called statically determinate when the number of interface degrees of freedom 
is equal to the ruraber of rigid bouy degrees of freedom of the structure at 
hand. Otherwise, the interface is called statically indeterminate. To show 
that ^ (20) is zero for a statically determinate 

interface, let us first state that the numerical values of the elements of 
this matrix are independent of the dynamical state of the structure. More 
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specifically, because the stiffness matrix is the same for P at rest 

or in motion, we can assume that the structure is in a state of static 
equilibrium without changing [ Therefore, without los 

of generality we consider the equilibrium equation of P under the action of 


KhW- 


N {^^} = ! 


\ I 

, using the partitioned form of I I Eq. (7), we can write 

[4n] j^i [*^i] j - |« : 

K] =!<1 


From Eq. (22) we can solve for \ P 

j ’'n 

\ p/ r..p r‘ r.,p 


K»r [<.] Hi - W 


where we used Eq. (13). Substituting Eq. (24) into Eq. (23) yields^ 


Kn] h* K]))-i| * |“t! 


At this point we should note that when the interface is statically 
determinate no stresses can be set up in P by the interface displacements 
I I'lndeed, for a statically determinate interface the matrix 
becomes a rigid body transformation, transforming the interface displacements 
into equivalent rigid body displacemnts of the non-intcrface degrees of 
freedom of P. Because, in addition we assumed tliat no t)tln‘r forces are acting 
on P, it is clear that | | is zero in Eq. (25), from which it follows tluit 
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This completes the proof. 


Finally, we note that in most applications of interest to this contract, 


the externally applied forces 
example, STS payloads will be eiiciosed 


t P 

enclot 


in tn 


and {F 




are nonexistent. For 


e cargo bay and will not be 


exposed to external forces. Therefore, we write 


N 




0 


(27) 


Taking into account Eqs- (19) and (27) we can now write the final form of 
Eq. (17). 




Wf J 









' 

1 




1 1 

- 1 - 


1 

CQ Ml 

X 1 
1 


0 1 



r ^ 

1 

O 

M 

' ~p ' 

‘K 


0 1 


in which 


{"b}- 


B 

B 


(28) 


(4) 


represents the generalized displacement vector of the free booster H- the 




vertor^J xj^ y contains all non-interface displacements of the booster 

""l I -P 


^represents the interface degrees of freedom. Furthermore, the 
vec tor jrepresents the non-*interface displacements of the payload P with 
respect to the interface. The matrices [Mfl] and [”p] represent the mass 


matrices of the booster and the payload respectively and 
represent the stiffness matrices. The matrix ^ 

a statically det 

.B 


is a transformation 

matrix, characteristic for the payload P. In case of a statically determinate 


interface, 




represents a rigid body transformation. Finally, 


N 


is 
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the externally applied force vector acting on the booster H* The 


matrices 


H N - [^p 7 N M 


contribute only to 


the interface degrees of freedom as can be seen from Eq. (20) and from 


[vl ["■•] 


0 0 


P P n 

^IN V 


L^FN 


IT 


0 S. 


0 T 


0 0 


0 M. 


(29) 


wi th 


K] - W'lKJH Wh] ^ [<.]■ 

The matrix essentially represents the payload mass 

transferred to the interface. Similarly, [^^1 represents the 

payload stiffness transferred to the interface. Note that when the interface 
is statically determinate no stiffness is transferred (£q. (26)). When the 
interface is statically indeterminate there are what is commonly called 


I ^ 

induce rigid body displacements in the payload but also strain's- These 


''constraint modes" [3], i.e. the interface displacements | not only 


strains cause the triple product 
2 ero. 


[hY H ["^p] 


to be different from 


Till?* ronrliidon r )t i n M<>t' f f on imi I 1 » 4 ‘ ioum of ttutl ion of I In* 

hooMtnr /pay I oad HysliMii. Thin i.il will ho uHotl (o d«*i iv<* .si*vi*i.il fnitln>»ls 

of solution of the equations of motion (28). 

3. SOLUTION OF THE EQUATIONS OF MOTION IN THE DISCRETE TIME DOMAIN 

As stated in the Introduction, the objective of this study is to determine 
design loads for the payload structure. These design loads are then used to 
calculate stresses and strains that would exist in the structural elements 
that make up the payload. The stresses and strains allow the designer to 
determin the correct physical and geometric properties of the elements (mass, 
stiffness, lengths, cross sections, etc.) so that the structure does not fail 
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when subjected to the external forces 
written as 


K] 


"ei'Wh] j«P 


An element loads equation is 


(31) 


which 
P 


i . 

e is the load vector of an individual element e of the payload 
} r T 




is the geometric 


> stiffness matrix of the element, and 

compatibility transformation. The vector time dependent 

displacement vector of the payload satisfying Eq. (28). Consequently, in 
order to determine | | need to solve Eq. (28). 


The most straightforward approach to determine | ^ is to solve Eq. (28) as 
a set of simultaneous second order differential equations. There are several 
well established response routines that handle such problems (Runge-Kutta, 
Newmark-Ch'in^Beta Numerical Integration, etc.). This Irect approach has the 
advantage of simplicity and accuracy. The obvious drawback is the high 
computational cost due to the large number of degrees of freedom used to 
describe today's aerospace models* Furthermore, this method does not take 
advantage of the fact that often only ^mall changes are made in the payload. 
H.:wever, thi method is still useful in the context of this study because it 
provides with reliable results that can be used for conqjarison purposes 
with other methods to be discussed shortly. 


fi. OP/rnje ANAI.YSIS [h,7] 

In this section we shall discuss a technique commonly kno*/n as modi 1 
analysis. l>.is approach will lead us to an alternate solution method for Eq. 
(18) '*1 we shall show that it has some definite advantages over the direct 
sol Lion of the set of differential equations (28) as discussed in Section 3. 


We start the process with the homogeneous set of equations extracted from 
the. top row of Eq. (3), 





( 31 ) 
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Associated with this equation is an eigenvalue problem, 

(-“b["bJ ^ [s]){*b} ■ {“} 

where -he vector| represents an eigenvector (mode shape) and an 
eigenvalue (natural frequency)* The solution of this eigenvalue problem 
essentially produces a linear transformation marix [^.l (modal transformation 
matrix) iu wliich I'ach column represent h a mode shape of the booster li. Tin* 
main property of this modal transformation is that in the new normal 
coordinate system^ , the equations of motion (31) become uncoupled, i.e* 

if we apply the transformation, 


to Eq. (31), and premultiply by 


[»b] 






we obtain 


where the coefficient matrices of 


W^WNKH'’ 


(33) 


(34) 


and ) q 




are now diagonal 


where 




written as 


[♦„]"["«][♦»] ■ [♦»]' [%][♦»] - [-b] 

I was normalized with respect to Pm 1 • Equation (3 

J [ Bj 


{y>) 

(34) can then be 


(36) 


The obvious advantage of applying the modal transformation Eq. (33) is 
that Eq. (36) now represents a set of decoupled independent second order 
differential equations that arc easily solved. The price to pay however, is 
the solution of eigenvalue problem (32). There are many well established 
eigenvalue problem ’'solvers** available (Jacobi, Rayleigh-Ri tz, etc.) 


The next step is to consider the homogenous equation, 

[M' N 


Note that, 


[^p] [^] [^p] “[”n’n]’ [^p]' [h] [s] = [Sn] 


( J7) 


(38) 
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So that Eq. (37) becomes 

T 


Klt^l ^ [q 


~P 


0 


(39) 


In the same way as we did for Eq. (31) we can introduce a modal 
transformation. 


with 


1 ■ [*U }‘'n 

K]K] = [']• [<^[<M] 


(40) 


(41) 


where we wrote instead of [ 1 to simplify the notation. The modal 

matrix columns the ”canti levered” mode shapes of the payload P, 

and J natural frequencies squared of the cantilevered payload 

(i.e. fixed interface) on its diagonal. Using Eqs. (40-41), we can write Eq. 
(39) as 


%l C“p] )%j ■ H 


(42) 


Let us now apply the following transformation to Eq. (28), 
B 


^ 1 





B f 1 

Jl__ ^ ' 



4> 1 

b1 

—— T“ “ 



I 

-P 

1 


! jP 
-T 


-p 


and premultiply by 
and (41) we obtain, 






(43) 


Taking into account Eqs. (35) 




Jb_ 


[“b] * ♦KVp*b| " 


i-!b_ 

-P^ T 1 

- ^ 1 ^ 

1 

-p 
‘«N 1 


" h]j 

1 

I 

p 

1 ''n / 


1 

I-'' 



N 

1 

1 

-!sJ 

0 

i f (A4) 

0 
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Note that | | can be partitioned as follows 

B 




(45) 


Equation (44) as it is, probably does not yield, any advantage over a 
Hrect solution of Eq. (28). However, in most practical applications there is 
a possibility of defining a so-called "cut-off frequency". In these cases a 


N 


shows that the energy 


Fourier series expansion of the force vector 

content of the high frequency components is smal I compared to the energy 
contained in the low frequency components. Practically, this means that the 

' I 

response of the structure due to the high frequency content of > can 

often be neglected. In this connection it should be noted that it is 


relatively difficult to excite the higher modes of the structure to any large 
extent, especially when | | contains a few elements (i.e. only a few 

application points). The idea then is to only retain these modes in Eq. (44) 
that have a frequency smaller than the cut-off frequency. This in turn, 
reduces the size of Eq. (44) considerably. Experience has shown that the 
introduction of a cut-off frequency is a workable concept. 


In conclusion, we can say that the introduction of a cut-off frequency 
leads to a reduction of the size of Eq. (44). Nevertheless, the solution of 
Eq. (44) for the modal displacement vector is still costly and again does not 
accommodate the special circumstances of small changes in the payload. A more 
serious problem however, is the representation of the interface in a model 
where a cut-off frequency is used. Indeed, in many cases a significant 
portion of the response is lost in the high frequency range due to the fact 
that the effect of the payload on the booster has a significant high frequency 
content. As will be shown in next section a need arises to improve the 
interface representation in the model for the booster. In sections 5 and 6 we 
shall discuss two methods that accomplish this. 
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5. THE RESIDUAL MASS AND STIFFNESS METHOD [S, 9, 10, 11-50, 5l] 


As explained in Section 4. in many cases it is possible to define a 
cut-off frequency which enables us to truncate the high modes in Eq. (44) 
thereby reducing the size of this equation* Obviously, some accuracy in the 
response of the structure is lost due to the truncation of these high modes. 
This loss of accuracy is especially apparent at the interface as we shall 
explain shortly. The residual mass and stiffness method, instead of omitting 
these high modes will replace them with a set of "residual modes". ‘Pin* 
computation of these residual modes does not require any knowledge of the 
payload so that they represent a one-time computation effort not to be 
repeated as long as the booster stays the same. In order to determine the 
residual inodes let us consider Eq. (33) 



(46) 


which represents the modal transformation for the booster B. Assuming a 
cut-off frequency was determined we can partition Eq. (46) as follows 




- 5 - 

H 


(47) 


where I T.. I represents the modes with frequencies less tfiaii l\ut cut-oil 
L “J r uT 


frequency 


and 

Hi 


those with higher frequencies. At this point one could 


r Hi L 

neglect I I and calculate the response as a linear combination of the lower 


mod 
loads J 


cs r^t^l only* Usually this yields a poor accuracy in the response and the 
dft L The reason 


The reason is that in most practical cases a significant part of the 
interface response is produced by the higher modes. Indeed, a typical 
interface is rather stiff and has little mass, i.e. that locally the interface 
has a high frequency content so that it responds significantly in the high 
frequency range. In truncating the higher modes the model does not include an 
adequate representation of that interface. The residual mass and stiffness 
method now, proposes Lo ret-ifn tlu‘ slat re eontr ihut ion to tlu‘ response* of I lu)^io 
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high modes* This leads to a much better representation of the interface. The 
static contribution can be obtained from the following static equation 


[•^] 


(48) 


derived from Eq. (3). Substituting Eq. (46) into Eq. (48) premultiplying 
recalling Eq. (35), yields 

] I’bI ■ 




(49) 


Because we are only interested in the high frequency part, let us write 
Eq. (49) as 

,T 




1 

0 1 
1 

L“b gj 






(50) 


So that from the bottom row in Ea. (50) we have 

IrB^ 

IT 


(51) 


[“b1 {<} - [♦“b1 

Finally, let us partitionj^ifg Jin non-interface and interface partitions, 




^BN 

’^BI 


(52) 


Substituting Eq. (52) into Eq. (51) we obtain, 

[<'] h“|- Kf Kzf 

In principal we can use Eq. (51) as it is and solve for 

,.-2-, r „ -.T I F® 

N 


(53) 




’S = [<] r 


(54) 
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which can be substituted in Eq» (47), yielding 




( 55 ) 


However, it should be noted that for every force component we keep, we add 

a degree-cf “freedom to the problem. If for example, » j contains many 

elements (i.e. many points of application) it may not pay off to use Eq.(55), 

i.e, we may as well keep all the modes in Eq. (46). If however | | 

contains a small number of elements (for example, in case of i landing or a 

rocket motor ignition) we can use Eq. (55) as it is, and obtain a much better 

response for few added degrees of freedom. However, because the cut-ot f ^ 

frequency was defined in such a way that all significant frequencies of 

are contained in the lower frequency range L, we can state that the booster 

model will adequately respond to \ i and no significant portion of the 

\ ^ / .(B). 

responni‘ will be lost. Therefore, we i «n omi I the term in <1*^ > in Eq . (51) 


The latter part 


altogether and jual keep the part in < l< . / 

1 B \ . ( B ) . ( ' ' 

in > ii* important because / will usually have a significant high 

frequency content (after all | | represents the effect of the payload on 

the booster and as such contains a wide range of frequencies). Because the 

interface usually has a high frequency content (as explained before) | | 

will induce a response at the interface primarily in the high frequency range 

which in turn will be, transmitted to the rest of the booster. 

On the other hand if i N i contains reaction elements due to some 
external contstraints (e.g. a dock) we wish to retain these elements as well 
because they are equivalent to elements of | j | sense that they 

represent the unknown effects of the constraints and also, the interface 
between the constraint and the booster usually has a high frequency content 
(e.g. connections between booster and dock). 
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Ignoring the term in | | i^i Eq. (53) and solving for yield* 


/ Hi 


r H"n 

r .H f 




L“b J 


(■'w 


( 56 ) 


which can be substituted into Bq. (47), yielding 

L 


{-b} ■[♦b| [♦“][“? JK] 




-2 

|v n 1 

r 1 

The term 
r Hi 

[♦«j 

r« J 

r»'J 


(57) 


ro placi 


represents the residual modes and they 
Also, note that modes only involve hoc ter quant itic^s 


which makes it a one-time computational (effort. 

Let us now derive the modally coupled equations of motion for the 
booster/payload system. First, we substitute Eq. (57) into the top row of Eq. 
(3) and then we premultiply by j"^L i , yielding 


» i 

f 


0 


L » i blM ] [*«] 


' 'D 4 l 


"L 

% 




][C] 




residual mass 


M 

1 

.[4] 

lf“f] 

[*"] 

I'J 


residual flexibility 




(58) 


Before proceeding, let us consider the homogeneous equation extracted from 
the lower half of Eq. (58) 

-4^ r V, (..B 


and solve the following eigenvalue problem 
leading to the modal transformation, 


(59) 


(60) 


(61) 
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and the properties 
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h]' W h]- M W 

We shall make use of these properties in deriving the modally couple! 
equations of motion of the booster/payload system. To this end let us write 
Eq. (28), 


L W? I -J 


.15. 


Kb + TpKpTp 


0 


T 

IpKpIp .. 


S, f 


-p 

X 


1 

N } 


) 


(03) 


Let us now introduce the following notations 




, r* 1 - [♦'•I [♦“ 1 1’^ r*“,f il 

I^BJ j KBII B I I BiJ R 


so that, combining Eqs. (57), (61) and (64) we can write. 


{"b} ■ [♦njl’B} 


We now define a transformation similar to Eq. (43) 


-p 

X 

n 


4> ! 

„^J 

I 

! -rP 




‘13 


where this time 


I ^ and (^g are given by Eq. (64) 


(64 ) 


(65) 


( 66 ) 
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The residual mass and Ntiffness technique essentially improves the 
interface representation in the model where the model is subjected to a 
frequency cut-off. 


The result Eq.(57) can then he used in any type of modal synthesis technique 
such as result Bq- (70). Due to this interface improvement it is now poHiiible 
to truncate the high hoostor modes while still obtaining an acceptable 
accuracy. As we shall discuss in section 9 9 the residual mass and si illness 
method •‘urns out to be the most efficient full-scale method currently 
aval latte in the literature. 

6 • niK MASS AND Si iyKNKSS I .OADlNt: TKCHN tgUE [ “i 2 , 1 0 , *> »J 


Another way of improving the interface representat ion in the booster mcuiel 
subject to frequency cut-ofl is given by a technique developed by Hrvida and 
Benfield and is based on Eq. (28) which we repeat here tor convenience, 


"b * 




1 I 


‘rV 




• j, 

"n 


^ ^ ■'pVp 




1 


1 ^ 1 

! 0 


1 0 1 



1 ^ J f 

n 



1 T. 

i 

P 1 1 { 




1 p p p 

1 

N 1 I I 



' 1 / 


(71) 


of solving eig«*nv«lue prohlein (32) Hruda and B**ufic>ld piopos-' to 
aolve the following eigenvalue problem. 


[«B*frv.'] ' [s ' • j"[ 


(72) 


yielding the modal trans format ion 


{s\ j-i; 


(7 i) 
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Again, the modal transformations (40) and (73) can be combined in 


B 



Substituting Eq. (74) into Eq. (71) and premul tip lying by 
yields. 


1 “ 
I 


0 ! 


(74) 

1 T 


1 T T p 

I j 


1 




_ I I 

1 

i 

+ 


L“pjJ 



(75) 


where we used the properties 







(76) 


Equation (75) now replaces Eq. (44). The main difference lies in the fact 
that in solving eigenvabie proWem (W), the booster interface is mass and 
stiffness loaded by I I I^P^^pj booster 

interface is loaded witn approximate dynamic effects of the payload. In doing 
so, the new modes frequencies will include a good 

representation of the interface. This allows us to reduce the number of 
booster modes in Eq. (75) according to the predetermined cut-off frequency. 

The disadvantage of this method in connection with the present study is chat 
eigenvalue problem (72) is dependent on the payload. This means that for 
every change in the payload we must solve this eigenvalue problem again 
although the booster does not change. This makes the Hruda/Benf i eld technique 
less suitable for our purposes. However, if the changes in P are smalJ, we 
can use the old booster modes as a first estimate to calculate th** new booster 
modes in a Raleigh-Ritz type eigenvalue problem solver. 


25 



OWGifJAL FAuc 
OF POOR OUALJTY 


7. THE COUPLED BASE MOTION TECHMIQUE [54.55,56,571 


The coupled base motion technique as presented in this section does not 
yield any immediate advantages over the methods presented in previous 
sections. However, it can be used as a starting point for possible short-cut 
methods (These possibilities will be investigated in Chapter II). In 
addition, this section will give us the opportunity to develop an alternative 
set of equations for Eq. (4). Indeed, we shall not only use "cantilevered'* 
displacements for the payload P but also for the booster B, while only the 
interface will be free. The derivation is very similar to the one in section 
2. Let us define a transformation similar to £q. (12) but now for the booster 


B, 


with 


^1 ■ H 

N " 'K] W 


(77) 


(78) 


Equation (14) can now be replaced by 




B 

1 


r !•’ 
1 
1 

1 


-B 

"n j 



B 

- 

0 ! 
1 

— 

I lo 


B 



1 p 

^ ' 

[ T% i 

1 

0 ! 
1 
1 

1 


T 

— P f 



\ P ' 

*i 


^ 0 I 

I !o _ 




Again, this transformation will eliminate the redundant set of 
displacements (^0) and in the process it will also eliminate the 

unknown reations' |rj | and | 

Introducing the notations 


- 



I ! ! 0 


I„ 1 T I 0 

1 B I 


B 1 B 1 

0 ! I I 0 

t 1 


1 ! 0 

1 1 

t 1 


1 1 
1 1 

0| Spj I 


0 ! ! 
1 1 

i 

O 

O 

1 


1 

o 


( 80 ) 
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where this time, 


Note that j different dimensions. 


have different dimensions. 


Substituting transformation (79) into Eq. (10) and premultiplying by A 


yields. 




^ Wp 


I i;Vp 


T T 
I 1 
FTP 


•>b'b 






ilv, ' ^ 


Equation (82) replaces Equation (28). 

The basic idea for a base drive method is the seperation of the booster 
response into two sepaiate parts 
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B 


= part due to the action 


of ^ 


(84) 


and 




HI 


- part due to the presence of the pc’ylo.'^d (^feedback) 


vdo) 


It is clear that vector (82) satisfies 




Wb Wb 


- 

"B 

1 

K 

+ 

‘X‘b ° 

1 

-B 

= * 

)‘P 

B B 

i 

“B 


_ » ^bVb. 

j 

B 

, 1 

T 

J 1 



1 

"i! 

! 1 

^b'^b 


( 86 ) 


The solution of Eq. (86) is a one-time computional effort because it only 
involves booster quantities. If we now substitute Eq. (83) into Eq. (82) and 
take into account Eq. (86) we obtain the following new set of equations 


‘X‘b : '^Vb . 

0 

i 

"n 


r^Xs; <> ! » 1 

( 

T ‘ T T , 

T 

) *BR 1 


— 

„ ' T T 1 

) brJ 

^b”b^b ’ ■'bVb^Wp 1 


) 

+ 

“ 1Tb*'b\%Vp' » 


■■ ■ “■ ^ — r 

1 T 1 


1 ^ 


n i n '.T 

f. T> 1 

0 -IpMpTp . 

IpMpIp J 

\ ^ i 


L 0 : 0 ,ipKpi^ 



( 0 j 1 

L-X. J 

1 jl.. * BF r 

T BF r 

j ^p^^p^i 1 "1 

Wp^'i 

J 

1 „ ‘ BF I 


( J 1 

( « I 


( 87 ) 
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where jxj^l and|xj*^| are known from Eq. 


BF 


( 86 ) 


The second idea of a base drive method is to consider the bottom partition 
of Eq. (87) and write it in the following form. 


[Ip«p‘p]t\l + * l?l) 


( 88 ) 


If one is only interested in the design of the payload^ Eqs. (88), (12) and 


/ **BR 

(83) is all we need, to deterains the response of the payload. | 

known we can %asc drive** the payload by the terms on the right-hand side of 

Eq. (88) to obtain • Of course, is coupled into the booster 

equations in Eq. (87). 


As mentioned before Eq. (87) does not yield any immediate advantages, but 

as will be discussed in the next chapter, Eq. (88) becomes very useful 
'BR 

i f { I (a feedback) is smal 1 . 


Equation (87) can also be %n:itten in terms of normal coordinates. To this 
end let us introduce the following transformation 



b 

where [<>n] «nd [ 4 .J are obtained from solving the following 
eigenvalue problems, 

[-‘“Bt'X'B] *['BVB]]t*Nl ■ i"l 
[-“? C^K‘B*’’JVp] ^ [’■bVb+^JVp]])*?! ■ 

with 



(89) 


(90) 

(91) 


( 92 ) 
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T 

■ B r^T. 


[-fi] [VBVWp]t<^i] = [I]. [<!>?] [tXv^pVp] ['*’?] 


(93) 


where again we used Che simpler notations instead of 


K3 


-B 1 and « B 




Substituting transformation (89) into Eq. (87) and premultiplying by the 
transpose of the square transformation matrix in Eq. (89) yields the modal 
form, 


I 


T 

rB^T. 






0 


bI 




1 T 

^P^Vn 




T / -BR 1 

■"bF 


0 «^nWp*i : 


(94) 
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8. A LOAD TRANSFORMATION CONSISTENT WITH MODAL SYCTHESIS TECHNIQUES [53] 


Before presenting the assessment section it will prove productive to 
discuss the topic of determining the internal structural loads in the payload 
members. As discussed in section 3, the reason for solving the equations of 
motion of the booster/payload system is the determination of the displacement 
vectorjxpl 3o that we can substitute this vector into Eq. (31), 

in order to obtain the internal structural loads | | on an individual member 

e of the payload P. In principal Eq. (95) could be used as it is, but this 
’’displacement** approach turns out to be very sensitive to inaccuracies in 
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|xp|; e*g. truncating high frequency modes as we did in section 5 could very 
easily lead to erratic loads Heuristically speaking, c ntains 

three parts, the static displacement, the rigid body displacement and the 
"vibrational" displacement. 


Therefore, if one has an error injxpjone necessarily affect the accuracy 
of all three parts. For this reason one prefers the so called "acceleration 
method". Basically this approach is capable of seperating the static and 
rigid bod> ^iarts from the "vibrational" displacement. As a consequence one 
only makes errors in l e "vibrational" part which often is the smallest part 
of the displaceirent vectorjXp Such an "acceleration" approach which is 
consistent with modal synthesis techniques was developed by Hruda and 
Jones. ^3 J 


Recalling Eq. (79) \.e can write 

Ivl - [fp] 


(96) 


so that from Eq. (95) we obtain 


(97) 


From the bottom roW of Eq. (82) we obtain 




and from the second row of eq. (82) we obtain 

’■?! ■ - [■'bVbK^I 




(98) 


(99) 


Expressions (98) and (99) can now be substituted into Eq. (97) yielding an 

p 

equation for in terms of accelerations. Many of the matrix 

multiplications involved in Eqs. (97-99) can be simplified by using a unit 
load solution, which is a feature of most finite element programs [53] 
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9- ASSESSMENT 


The methods as discussed in sections 1-8 are believed to be the currently 
most prominent full-scale methods* Some of these are improvements or 
adaptations of previously existing approaches (Hurty, MacNeal, Bamford, 
Craig/Bampton, etc*) A study by R. Hruda [il-50] showed that the residual 
mass and flexibility approach is the most effective in terms of cost and 
convenience* As a test structure, Hruda used two planar trusses cc “ ed 
together by a statically indeterminate interface (Figure 2.) Five different 
techniques were compared to the exact solution, i.e. the solution in the 
discrete time domain as discussed in section 3*: 

1* Hruda/Benfield Technique (section 6): inertial coupling of truss-2 
constrained modes onto free-free modes of truss-1 which was mass and 
stiffness loaded at its truss-2 interface degrees-of-freedom by the 
interface properties of truss-2. (IMSL) 

2. Craig/Bampton Technique (modal version of Eq. (82)): inertial 
coupling of truss-1 and truss-2 constrained modes onto a free-free 
modal representation of the interface degrees-of-f reedom. ( I/F) 

3. MacNeal Technique : residual flexibility approach of coupling truss-2 
constrained modes onto free-free modes of truss-1 which creates 
stiffness coupling (residual mass not included) . (RFSWOM) 

Rubin Technique (the residual mass and flexibility technique - 
section 5): coupling of truss-2 constrained modes onto free-free 
modes of truss- 1 which yields only inertial coupling, and, by 
consistent application to the mass and stiffness terms in the 
equations of motion, yields both residual stiffness and residual mass 
terms. (RPIWM) 
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5. Rubin Technique but without residual nass contribution for truss-l* 

Tne truss problem as illustrated in Figure 2. represents a 
planar problem with three rigid body degrees of freedom (two 
translational 9 one rotational). Each pinned joint has two 
translational degrees of freedom. The interface is statically 
indeterminate because there are six interface degrees of freedom. 

The heavy masses (asymmetric with respect to interface) are added to 
produce interface distortion. The forcing function is a ramp 
function. (RFIWOM) 

The **ex«act” results, against which all comparisons were made, 
were obtained by extracting eigenvalues, eigenvectors, and loads 
directly from a finite element discrete/physical model using no modes 
at all. 

Five different cases we investigated 
EXACT: Discrete modal 70 DOF 

CASE A: Modally coupled, 70 modes retained (*10GZ) 

CASE B: Modally , coupled, 50 modes retained (=71Z) 

CASE C: Modally coupled, 19 modes retained (*46Z) 

CASE F: Modally coupled, 19 modes retained (=27%) 

Hruda used the following comparison values: 

Frequencies: percent error against the **exact” solution. 
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Modes: 


An error vector is formed ) end its norm is 

calculated (which is defined as the Root Square Sum of 
the elements of the vector); the comparison value is 
then defined as the norm of the error divided by the 
norm of the base/exact mode. Note that the norms are 
based on the modal amplitudes of all degrees of 
freedom from the coupled system. 

Loads: Loads were calculated at the truss interface on both 

the trus8“l and truss-2 joints. A percent error of 
the absolute value of the largest (either raaximuu or 
minimum) value from a given case against the absolute 
v.iluo of tl»c lar^i'Ht valm* from lln' c'xact soltiCion. 


i.e. 


Frequency comparison value » 



X 100 


Mode comparison value 
Load comparison value 


RSS ( (|> 


N 


RSS ^ 


E 


X 100 


X 100 


where E^Exacti and N«Case being compared. 

The results are presented in Tables 1-12 
For the 100% case-Ay the MacNeal technique requires the inversion of the 
residual flexibility matrix to obtain a ’’residual stiffness”. When 
attempting to retain all (100%) of the modes* this residual flexibility 
matrix is a function of the interface highest frequency modal amplitudes 
which can cause an ill-conditioned matrix (as in the present case). Since 
this is an unrepresentative case* it should not be deduced that this is an 
unacceptable technique. As can be seen in succeeding cases* where more 
residual modes are available* the MacNeal technique falls into line with 
other techniques. Note that in cases B,C* and F, in both the frequency 
and mode shape comparisons* that the MacNeal and the Rubin technique 
without residual mass are identical* thereby numerically supporting R. 
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Coppolino's contention that these two techniques are equivalent for modal 
synthesis. In comparing the loads it is seen that the MacNeal column for 
case~A reflects the propagation of the ill-conditioning mentioned 
earlier. Loads were calculated by the modal acceleration technique 
(section 8) for all methods except for the MacNeal technique. 

Due to the stiffness coupling involved in the MacNeal method, a complete 
modal acceleration technique for calculating loads could not be used, 
therefore, the modal displacement techniques of calculating loads was used. 
Because of this, the larger loads inaccuracies for this method must be 
attributed to the method of loads calculation and not to the method it^ielf. 

In conclusion we can state that methods I through 5 are acceptable. 
However, the Rubin Technique (Residual Mass and Stiffness Approach) seems to 
outweigh the other approaches in terms of cost and convenience. Again, it 
should be noted that this method does not require any knowledge of payload 
properties which makes it very valuable for analysis of STS-applications. 



Figure 2 Structure Used for Comparing Coupling Techniques 
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Table 1 : Frequency Comparison 

Case A » 100% of Available Modes 
% Dift' for Various Modal Coupling Tech 


Mode 

ho 

Exact Freq 
(Hz) IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

.76 

-.00 

-.00 

-49.17 

-.00 

-.00 

5 

1.75 

-.00 

-.00 

6.75 

-.00 

.00 

6 

2.84 

-.00 

-.00 

-21.81 

-.00 

-.00 

7 

3.08 

-.00 

-.00 

-3.63 

-.00 

.00 

8 

3.80 

-.00 

-.00 

-9 . 24 

-.00 

.00 

9 

4.62 

-.00 

-.00 

-16.35 

-.00 

.00 

10 

5.11 

-.00 

-.00 

-24.08 

-.00 

.00 

11 

5.50 

-.00 

-.00 

-6.02 

-.00 

.00 

12 

5.81 

-.00 

-.00 

-3.31 

-.00 

.00 

13 

7.69 

-.00 

-.00 

-6.08 

-.00 

.00 

14 

8.69 

-.00 

-.00 

-.60 

-.00 

.00 

15 

9.14 

-.00 

-.00 

-3.32 

-.00 

.00 

16 

9.42 

-.00 

-.00 

-1.76 

-.00 

.00 

17 

5.73 

-.00 

-.00 

-.92 

-.00 

.00 

18 

9.85 

-.00 

-.00 

-.98 

-.00 

-.00 

19 

10.36 

-.00 

-.00 

-1.36 

-.00 

.00 

20 

10.43 

-.00 

-.00 

-.33 

-.00 

.00 

21 

10.79 

-.00 

-.00 

-2.99 

-.00 

.00 

22 

10.90 

-.00 

-.00 

-2.55 

-.00 

.00 

23 

11.37 

-.00 

-.00 

-1.89 

-.00 

.00 

24 

11.49 

-.00 

-.00 

.80 

-.00 

.00 

25 

11.78 

-.00 

-.00 

.39 

-.00 

.00 

26 

11.96 

-.00 

-.00 

-.12 

-.00 

.00 

27 

12.03 

-.00 

-.00 

.19 

-.00 

.00 

28 

12.20 

-.00 

-.00 

.40 

-.00 

.00 

29 

12.43 

-.00 

-.00 

.12 

-.00 

.Oo 

30 

12.50 

-.00 

-.00 

5.58 

-.00 

.00 

31 

12.75 

-.00 

-.00 

4.00 

-.00 

.00 

32 

13.29 

-.00 

-.00 

2.54 

-.00 

.00 

33 

13.51 

-.00 

-.00 

3.12 

-.00 

.00 

34 

14.22 

-.00 

-.00 

-.19 

-.00 

.00 

33 

14.53 

-.00 

-.00 

.31 

-.00 

.00 

36 

14.86 

-.00 

-.00 

-.11 

-.00 

.00 

37 

15.19 

-.00 

-.00 

1.59 

-.00 

.00 

38 

15.54 

-.00 

-.00 

1.23 

-.00 

.00 

39 

15.69 

-.00 

-.00 

2 50 

-.00 

.00 

40 

16.16 

-.00 

-.00 

.55 

-.00 

.00 

41 

16.17 

-.00 

-.00 

.66 

-.00 

.00 

42 

16.28 

-.00 

-.00 

3.34 

-.00 

,00 

43 

16.86 

-.00 

-.00 

.46 

-.00 

.00 

44 

17.07 

-.00 

-.00 

.21 

-.00 

.00 

45 

17.17 

-.00 

-.00 

5.12 

-.00 

.00 

46 

18.17 

-.00 

-.00 

,03 

-.00 

.00 

47 

18.17 

-.00 

-.00 

10,37 

-.00 

.00 
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Table 1: Frequency Ccnpar 1 son * (Cone 1) 

Case A « lOOZ Available Modes 
Z Diff for Varous Modal Coupling Tech 


Mode 

No 

E’^act Freq 
(Hz) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

48 

20.11 

-.00 

-.00 

.81 

-.00 

.00 

49 

20.27 

-.00 

-.00 

3.02 

-.00 


50 

21.11 

-.00 

-•00 

.14 

-.00 

.00 

51 

21.15 

-.00 

-.00 

.08 

-.00 

.00 

52 

21.26 

-.00 

-.00 

.07 

-.00 

.00 

53 

21.29 

-.00 

-.00 

.53 

-.00 

.00 

54 

21.45 

-.00 

-.00 

.02 

-.00 

.00 

55 

21.47 

-.00 

-.00 

2.91 

-.00 

.JO 

56 

22.06 

-.00 

-.00 

.45 

-.^;o 

.00 

57 

22.15 

-.00 

-.00 

2.18 

-.00 

.00 

58 

22.53 

-.00 

-.00 

5.09 

-.00 

.00 
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Table 2: Frequency Comparison 

Case B 3 71Z of Available Modes 
X Diff for Various Modal Coupling Tech 


Mode 

No 

Exact Freq 
(Ha) 

IMSL 

I/F 

RFSWOM 

F’ IWM 

RFIWOM 

4 

.76 

.00 

.00 

.00 

.00 

.0^ 

5 

1.75 

.00 

.00 

.00 

.00 

0 

6 

2.84 

.00 

.00 

.00 

.00 

02 

7 

3.08 

.00 

.00 

.00 

.00 

.06 

8 

3.80 

.00 

.00 

.00 

.00 

.00 

9 

4.62 

.00 

.00 

.04 

.00 

.04 

?0 

5.11 

.00 

.00 

.00 

.00 

.00 

11 

5.50 

.00 

.00 

.05 

.00 

.05 

12 

5.81 

.04 

.05 

.03 

.03 

.03 

13 

7.60 

.00 

.00 

.21 

.00 

.21 

14 

8.69 

.02 

.02 

.03 

,01 

,03 

15 

9.14 

.04 

.06 

.12 

.04 

.12 

16 

9.42 

.01 

.01 

.01 

.00 

.01 

17 

9.73 

.01 

.02 

.02 

.01 

.02 

18 

9.85 

.03 

.04 

.03 

.03 

.03 

19 

10.36 

.01 

.01 

.06 

.01 

.06 

?0 

10.43 

.00 

.00 

.01 

.00 

.01 

2i 

10.79 

.03 

.05 

.04 

.02 

.04 

33 

10.90 

.03 

.03 

.04 

.02 

.04 

23 

11.37 

.00 

.00 

.03 

,00 

.03 

24 

11.49 

.00 

.01 

.05 

.00 

.05 

25 

11.78 

.03 

.03 

.08 

.03 

.08 

26 

11.96 

.08 

.09 

.08 

.07 

.08 

27 

12.03 

.02 

.03 

.04 

.02 

.04 

28 

12.20 

.13 

.13 

.10 

.10 

.10 

29 

12.43 

.01 

.01 

.02 

.00 

.02 

30 

12.50 

.03 

.09 

.07 

.03 

.07 

31 

12.75 

.07 

.07 

.08 

.06 

.08 

32 

13.29 

.04 

.05 

.06 

.02 

.06 

33 

13.51 

.03 

.07 

.07 

.05 

.07 

34 

14.22 

.13 

.13 

.78 

.09 

.78 

35 

14.53 

.09 

.13 

.29 

.07 

.29 

36 

14.86 

.01 

.01 

.17 

.01 

.17 

37 

15.19 

.06 

.07 

.11 

.05 

.11 

38 

•5.54 

.01 

.03 

1.17 

.04 

1.17 

39 

13.69 

.02 

.02 

1.61 

.03 

1.61 

40 

16.16 

.01 

.02 

.13 

.09 

.13 

41 

16.17 

.07 

.08 

2.68 

.12 

2.68 

42 

16-28 

.37 

.44 

3.88 

.34 

3.38 

43 

16.86 

.13. 

.12 

7.73 

.25 

7.73 

44 

17.07 

.06 

. 10 

18.51 

1 .52 

18.51 

45 

17.17 

1 .00 

1 .00 

60.94 

5 . 86 

66.94 

46 

18.17 

.02 

.02 

72.29 

9.54 

72.29 
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Table 2: Frequency Comparison (Concl) 

Case B ” 71Z of Available Modes 
X Diff for Various Modal Coupling Tech 


Mode 

No 

Exact Freq 
(Hz) 

IMSL 

1/F 

Rt'SUOM 

RFIWM 

RFIWOM 

47 

18.17 

.03 

.04 

122.22 

11.51 

122.22 

48 

20.11 

.07 

.09 

473.05 

7.55 

473.05 

49 

20.27 

.31 

.3 

962.70 

8.65 

902.70 

50 

21.11 

4.07 

4.20 

536.67 

19.40 

536.67 
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Table 3: Frequency Comparison 

Case C = 46X of Available Modes 



Z 

Diff for 

Various Modal 

Coupling Tech 



Mode 

No 

Exact Freq 
(Hz) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

,76 

• 00 

.00 

.02 

.00 

.02 

5 

1.75 

.02 

.00 

.00 

.00 

.00 

6 

2.84 

.04 

.02 

.12 

.01 

.12 

7 

3.08 

.01 

.00 

.28 

.00 

.28 

8 

3.80 

.02 

.00 

.00 

.00 

.00 

9 

4.62 

.12 

.05 

.26 

.02 

.28 

10 

5.1’ 

.02 

.01 

.01 

.00 

.01 

11 

5.50 

.26 

.10 

.35 

.05 

.35 

12 

5.81 

.16 

.24 

.12 

08 

. 12 

13 

7.69 

.03 

.04 

1.55 

.06 

1.55 

14 

8.69 

.14 

.22 

.17 

.08 

.17 

15 

9.14 

.32 

.33 

.41 

.17 

.41 

16 

9.42 

.05 

.07 

.07 

.01 

.07 

17 

9.73 

.13 

.18 

.34 

.08 

.^4 

IS 

9.85 

.30 

.27 

.31 

.18 

.31 

19 

10.36 

.11 

.11 

.51 

.12 

-51 

20 

10.43 

.09 

.04 

.45 

.02 

.45 

21 

10.79 

.17 

.33 

.73 

.11 

.73 

22 

10.90 

.41 

.42 

3.43 

.55 

3.43 

23 

11.37 

.06 

-04 

3.23 

1.49 

3.23 

24 

11.49 

.04 

.11 

5.14 

3.66 

5.14 

25 

11.78 

.17 

.17 

5.52 

2.76 

5.52 

26 

11.96 

.40 

.90 

4.79 

3.70 

4.79 

27 

12.03 

.30 

3.16 

53.65 

4.71 

53.65 

28 

12.20 

1.73 

2.35 

77.63 

3.91 

77.63 

29 

12.43 

.41 

4.66 

80.30 

14.58 

80.30 

30 

12.50 

1.94 

12.32 

143.04 

22.38 

143.04 

31 

12.75 

3.52 

14.90 

167.76 

29.91 

167.76 

32 

13.29 

11.88 

61.19 

630.91 

63.27 

630.19 
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Table 4: Frequency Comparisor 

Case F « 27Z of Available Modes 




X 

Diff for 

Various Modal Coupling Tech 



Mode 

No 

Exact Freq 
(Hx) 

IliSL 

l/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

.76 


.00 

.00 

.05 

.00 

.05 

5 

1.75 


.02 

.00 

.02 

.00 

.02 

6 

2.8A 


.06 

.06 

.31 

.03 

.31 

7 

3.08 


.07 

.DO 

4.32 

.19 

4.32 

8 

3.80 


.03 

.01 

.03 

.00 

.03 

9 

4.62 


.18 

.14 

.71 

.08 

.71 

10 

5.11 


.05 

.03 

.03 

.01 

.03 

11 

5.50 


.46 

.34 

1.64 

.49 

:.64 

12 

5.81 


:• .03 

1.13 

23.29 

4.70 

23.29 

13 

7.69 


.14 

.16 

15.39 

8.20 

15.39 

14 

8.69 


1.05 

1.14 

45.28 

1.30 

45.28 

15 

9.14 


1.20 

6.85 

84.48 

8.49 

84.48 

16 

9.42 


.37 

27.17 

98.00 

27.57 

98.00 

17 

9.73 


2.23 

36.67 

97.08 

36.41 

97.08 

18 

9.85 


18.25 

37.17 

203.74 

37.43 

203.74 

19 

10.36 


22.07 

103.62 

279.52 

105.23 

279.52 
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^able ^ : Mode Shape Comparison 

Case A * 1001 of Available Modes 



Z 

Diff for 

Various Modal 

Coupling Tech . 



Mode 

No 

Exact Freq 
( HZ ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.00 

.00 

104.09 

.00 

.00 

5 

13.95 

.00 

.00 

92.76 

.00 

.01 

6 

13.45 

.00 

.00 

104.91 

.00 

.01 

7 

13.70 

.00 

.00 

94.49 

.00 

.02 

8 

12.23 

.00 

.00 

132.95 

.00 

.02 

9 

11.21 

.00 

.00 

100.73 

.00 

.03 

10 

9.54 

.00 

.00 

153.35 

.00 

.05 

11 

8.75 

.00 

.00 

104.39 

.00 

.12 

12 

13.25 

.00 

.00 

129.05 

.00 

.08 

13 

10.51 

.00 

.00 

186.06 

.00 

.18 

14 

11.49 

.00 

.00 

96.78 

.00 

.04 

15 

11.24 

.00 

.00 

102.21 

.00 

.16 

16 

14.5 

.00 

.00 

71.19 

.00 

.07 

17 

13.16 

.00 

.00 

141.30 

.00 

.09 

18 

13.99 

.00 

.00 

90.01 

.00 

.02 

19 

15.47 

.00 

.00 

81.01 

.00 

.12 

20 

13.49 

.00 

.00 

99.38 

.00 

.05 

21 

12.74 

.00 

.00 

136.40 

.00 

.13 

22 

15.25 

.00 

.00 

138.19 

.00 

.15 

23 

16.56 

.00 

.00 

132.20 

.00 

.15 

24 

18.06 

.00 

.00 

100.18 

.00 

.25 

25 

17.70 

.00 

.00 

132.72 

.00 

.18 

26 

17.29 

.00 

.00 

142.58 

.00 

.17 

27 

16.94 

.00 

.00 

68.01 

.00 

.09 

28 

17.58 

.00 

.00 

69.89 

.00 

.03 

29 

17.06 

.00 

.00 

41.89 

.00 

.19 

30 

16.31 

.00 

.00 

138.22 

.00 

.28 

31 

15.86 

.00 

.00 

149*30 

.00 

.01 

32 

18.16 

.00 

.00 

138.85 

.00 

.08 

33 

16.92 

.00 

.00 

148.49 

.00 

.11 

34 

15.21 

.00 

.00 

91.84 

.00 

.84 

35 

18.67 

.00 

.00 

75.41 

.00 

.52 

36 

18.48 

.00 

.00 

49.75 

.00 

.36 

37 

17.23 

.00 

.00 

89.32 

.00 

.16 

38 

18.14 

.00 

.00 

145.84 

.00 

1.36 

39 

18.12 

.00 

.00 

146.47 

.00 

.72 

40 

18.82 

.00 

.00 

160.19 

.00 

.51 

41 

19.22 

.00 

.00 

126.88 

.00 

.56 

42 

J5.63 

.00 

.00 

162.42 

.00 

.90 

43 

18.47 

.00 

.00 

124.22 

.00 

1.12 

44 

19.44 

.00 

.00 

61.05 

.00 

.45 

45 

16.54 

.00 

.00 

135.41 

.00 

,84 

46 

18.17 

.00 

.00 

54.79 

.00 

.32 
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Table 5: Mode Shape Compariaon (Concl) 

Case A - lOOZ of Available Modes 
X Diff for Various Modal Coupling Tech. 


Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

47 

18.33 

.00 

.00 

142.35 

.00 

.30 

48 

18.25 

.00 

.00 

145.79 

.00 

2.22 

49 

19.16 

.00 

.00 

139.28 

.00 

.27 

50 

19.29 

.00 

.00 

65.342 

.00 

2.16 

51 

19.51 

.00 

.00 

62.72 

.00 

1.22 

52 

19.17 

.00 

.00 

52.42 

.00 

3.43 

53 

18.89 

.00 

.00 

105.50 

.00 

4.69 

54 

19.52 

.00 

.00 

38.21 

.00 

1.48 

55 

19.13 

.00 

.00 

103.39 

.00 

2.35 

56 

11.17 

.00 

.00 

161.11 

.00 

5.68 

57 

19.81 

.00 

.00 

11.29 

.00 

.17 

58 

19.81 

.00 

.00 

135.33 

.00 

.09 
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Table 6: Mode Shape Comparison 

Case B * 71Z of Available Modes 



J 

! Diff for 

Various Mddal 

Coupling Tech. 



Mode 

No 

Exact Freq 

(HZ) IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.00 

.00 

.00 

.00 

.00 

5 

13.05 

.02 

.01 

,.01 

.01 

.01 

6 

13.45 

.04 

.03 

.18 

.02 

.18 

7 

13.70 

.07 

.01 

.23 

.01 

.23 

8 

12,23 

.08 

.06 

.06 

.05 

.06 

9 

11.21 

.12 

.08 

. 33 

.Ob 

.33 

10 

9.54 

.16 

.18 

. 16 

.15 

.16 

11 

8.75 

.22 

.22 

.66 

.14 

.66 

12 

13.25 

.80 

.89 

.7o 

.70 

.76 

13 

10.51 

.30 

.41 

1.78 

.27 

1.78 

14 

11.49 

.91 

).21 

.88 

.57 

.88 

15 

11.24 

1.75 

2.05 

2.53 

1.67 

2.53 

16 

14.55 

.61 

.71 

.95 

.34 

.95 

17 

13.16 

1.02 

1.10 

1.02 

.70 

1.02 

18 

13.99 

1.52 

1.67 

1.63 

1.50 

1.63 

19 

15.47 

1.28 

1.03 

1.86 

.77 

1.86 

20 

13.49 

1.30 

1.03 

1.21 

.75 

1.21 

21 

12.74 

1.54 

2.02 

2.20 

1.23 

2.20 

22 

15.25 

1.93 

1.84 

1.98 

1.50 

1.98 

23 

16.56 

.38 

.62 

3.31 

.12 

3.31 

24 

18.06 

.32 

.81 

3.56 

.31 

3.56 

25 

17.70 

3.44 

3.63 

4.94 

3.43 

4.94 

26 

17.29 

9.71 

9.84 

8.9i 

8.36 

8.91 

27 

16.94 

8.37 

8.24 

7.18 

7.04 

7.18 

28 

17.58 

7.16 

7.15 

6.09 

5.94 

6.09 

29 

17.06 

1.11 

2.22 

3.50 

.79 

3.50 

30 

16.31 

3.83 

5.48 

5.80 

3.52 

5.80 

31 

15.88 

3.90 

3.89 

4.95 

3.86 

4.95 

32 

18.16 

2.81 

2.72 

4.45 

2.02 

4.45 

33 

16.92 

4.06 

3.29 

3.91 

2.92 

3.91 

34 

15.21 

6.36 

4.34 

25.34 

4.43 

25.34 

35 

18.67 

5.53 

5.12 

22.71 

4.42 

22.71 

36 

18.48 

1.64 

1.62 

16.67 

1.79 

16.67 

37 

17.23 

4.65 

5.05 

12.42 

3.98 

12.42 

38 

18.14 

1.43 

2.71 

71.43 

6.38 

71.43 

39 

18.12 

2.45 

2.06 

189.04 

4.32 

189.04 

40 

18.82 

35.48 

19.60 

104.85 

60.42 

104.85 

41 

19.22 

39.52 

27.36 

128.81 

60.10 

128.81 

42 

15.63 

26.78 

27.93 

126.75 

28.17 

126.75 

43 

18.47 

9.51 

9.90 

140.52 

38.29 

140.52 

44 

19.44 

14.36 

15.85 

144.04 

126.68 

144.04 

45 

16.54 

24.08 

24.65 

213.55 

142.78 

213.55 

46 

18.17 

3.25 

4.86 

150.22 

123.17 

150.22 
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Table 6: Mode Shape Comparison (Concl) 

Case B ■ 71Z of Available Modes 



X 

Diff for 

Various Modal 

Coupling Tech* 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

i\i 

18.33 

4.54 

5.77 

250.51 

144.13 

250.51 

48 

18.25 

62.46 

70.50 

553.15 

77.72 

553.15 

49 

19.16 

57.97 

66.55 

965.00 

130.07 

985.90 

50 

19.29 

114.23 

116.94 

930.68 

128.01 

930.68 
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Table 7: Mode Shape Comparison 

Case C * 46Z of Available Modes 
X Diff for Various Modal Coupling Tech, 


Mode Exact Freq 


No 

(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

• 02 

.02 

.03 

.02 

.03 

5 

13.95 

• 18 

.04 

.03 

.03 

.03 

6 

13.45 

.47 

.45 

.49 

.31 

.49 

7 

13.70 

• 3i 

.10 

.55 

.13 

.55 

8 

12.23 

.54 

.23 

.18 

.16 

.18 

9 

11.21 

• 1.59 

1.30 

2.00 

.90 

2.00 

10 

9.54 

1.10 

.74 

.41 

.36 

.41 

11 

8.75 

2.79 

3.01 

3.61 

1.79 

3.61 

12 

13.25 

2.50 

3.48 

2.42 

1.79 

2.42 

13 

10.51 

1.92 

2.35 

9.62 

2.83 

9.62 

14 

11.49 

4.84 

7.14 

5.19 

3.88 

5.19 

15 

11.24 

8.72 

9.37 

8.63 

5.92 

8.63 

16 

14.55 

3.33 

4.77 

3.71 

1.62 

3.71 

17 

13.16 

9.03 

90.46 

13.85 

5.90 

13.85 

18 

13.99 

11.17 

10.37 

13.28 

7.67 

13.28 

19 

15.47 

12.62 

U.13 

68.73 

9.12 

68.73 

20 

13.49 

13.23 

10.81 

31.34 

7.28 

81.34 

21 

12.74 

7.56 

16.01 

46.68 

10.40 

46.68 

22 

14.25 

16.60 

19.11 

54.62 

20.68 

54.62 

23 

16.56 

6.01 

5.57 

123.28 

97.95 

123.28 

24 

18.06 

4.87 

9.79 

129.15 

145.42 

129.15 

25 

17.70 

17.82 

16.45 

135.79 

155.65 

135.79 

26 

17.29 

68.71 

127.43 

150.09 

124.74 

150.09 

27 

16.94 

72.98 

144.11 

134.45 

120.45 

134.45 

28 

17.58 

103.74 

98.65 

152.44 

150.49 

152.44 

29 

17.06 

81.82 

139.33 

166.49 

128.91 

166.49 

30 

16.31 

110.88 

130.70 

145.30 

111.5S 

145.30 

31 

15.86 

95.31 

153.24 

243.62 

141.77 

243.62 

32 

18.16 

101.15 

122.19 

695.80 

128.54 

695.89 
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Table 8: Mode Shape Comparison 

Case F « 27Z of Available Modes 



1 

Diff for 

Various Modal 

Coupling Tech* 



Mode 

No 

Exact Freq 
(HZ) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

4 

14.84 

.04 

.05 

i06 

.03 

,06 

5 

13.95 

.24 

.07 

.26 

.06 

.26 

6 

13.45 

.71 

.89 

2.63 

.71 

2.63 

7 

13.70 

1.02 

.20 

8.69 

2.36 

8.69 

8 

12.23 

.74 

.39 

1.48 

.42 

1.48 

9 

11.21 

2.58 

2.75 

4.62 

1.96 

4.62 

10 

9.54 

2.84 

2.19 

1.40 

.98 

1.40 

11 

8.75 

10.05 

9.34 

38.19 

23.43 

38.19 

12 

13.25 

10.96 

11.05 

65.52 

28.35 

65.52 

13 

10.51 

7.03 

7.59 

128.49 

53.35 

128.49 

14 

11.49 

29.41 

33.79 

159.45 

43.11 

159,45 

15 

11.24 

33.59 

66.95 

167.22 

88.27 

167.22 

16 

14.55 

30.24 

113.08 

161.12 

119.44 

161.12 

17 

13.16 

111.59 

138.72 

161.33 

140.08 

161.33 

18 

13.99 

123.87 

130.93 

148.55 

132.90 

148.55 

19 

15.47 

148.78 

127.29 

379.09 

127.01 

379.09 
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Table 9 

Comparisons of Maximua' Absolute Values of Interface Loads 
Case A • lOOZ of Available Modes 


Load No. 

i^ixact Load (lbs) 

IMSL 

Percent 

I/F 

Difference - 
RFSWOM 

Abo. Max. 
RFIWM 

Loads 

RFIWOM 

1 

*481.999 

0 . 

-.00 

3419.56 

.00 

-.05 

3 

-202.138 

0. 

-.00 

10154.67 

.00 

-.08 

5 

-498.819 

0 . 

.00 

3957.31 

-.00 

.11 

7 

474.713 

0. 

-.00 

2208.71 

.00 

-.03 

9 

191.901 

0. 

-.00 

22733.86 

.00 

-.02 

11 

486.870 

0 . 

.00 

17322.92 

.00 

.03 


IMSL * Inertial Coupling W/ Mass and Stiffness Loading 
l/F = Interface Method of Inertial Coupling 

RFSWOM = Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIWM * Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM = Residual Flexibility with Inertial Coupling, without Residual Mass 


Table 10 

Comparisons of Maximum Absolute Values of Interface Loads 
Case B “ 71Z of Available Modes 


Percent Difference - ABS. Max. Loads 


Load 

No. Ex cc Load (lbs) 

IHSL 

I/F 

RFSWOM 

RFIWM 

RFIWOM 

1 

-481.999 

.10 

-.00 

2.08 

. 10 

-.38 

3 

-202.138 

.26 

.27 

2.62 

.27 

.26 

5 

-498.819 

.12 

.13 

-.13 

.14 

-.20 

7 

474.713 

.17 

.15 

-1.36 

.18 

.22 

9 

191.901 

.15 

.18 

.26 

.13 

-.24 

11 

486.870 

.19 

.19 

2.00 

.24 

-.10 

IMSL 

» Inertial Coupling W/ 

Mass and 

Stiffness 

Load i ng 




I/F = Interface Method of Inertial Coupling 

RFSWOM * Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIWM * Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM * Residual Flexibility with Inertial Coupling, without Residual Mass 
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Table 11 

Comp ax Isons of Maximum Absolute Values of Interf< Loads 
Case C ■ 46X of Available Modes 

Percent Differerce - ABS. Max. Loads 


Load No. 

Exact Load (lbs) 

IMSL 

I/F 

RFSWOM 

RFIWM 

RFIWC 

1 

-481.999 

-1.32 

-1.29 

-4.07 

-1 .09 

,29 

3 

-202.138 

2.34 

2,33 

3,10 

2.43 

3 . ! 9 

5 

-498.819 

.25 

,39 

-7.33 

.45 

.17 

7 

474.713 

-1.48 

-1,43 

2.93 

-1.36 

-.34 

9 

191 901 

.00 

-.12 

-25.44 

-.07 

-1.48 

11 

486.870 

-.85 

-.77 

-3,88 

-.78 

-.73 


IMSL * Inertial Coupling W/ Mass and Stiffness Loading 
I/F = Interface Method of Inertial Coupling 

RFSWOM = Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIWM = Residual Flexibility with Inertial Coupling, wiwh ^esidual Mass 
RFIWOM * Residual Flexibility with Inertial Coupling, without Residual Mass 
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Table 12 

Comparisons of Maximum Absolute Values of Interface 
Case F - 27% of Available Modes 

Percent Difference - ABS. Max. Loads 


Load No. 

Exact Load (lbs) 

:msl 

I/F 

rfswom 

RFIWM 

RFIWOM 

1 

-481.999 

-2.45 

-'•.88 

6.38 

-6.66 

-2.60 

3 

-202.138 

3.04 

-2.85 

5.68 

-10.99 

.27 

5 

-498.819 

-2.24 

-2.61 

8.02 

-4.49 

-1.85 

7 

474.713 

-3.82 

-4.44 

10.81 

-8.79 

-10.31 

9 

191.901 

-2.79 

-.54 

,21.52 

-7.54 

-14.78 

11 

486.870 

-.08 

-.48 

17. Of 

-4.28 

-7.60 


IMSL = Inertial Coupling W/ Macs and Stiffness Loading 
I/F = lntcr> *co Method of Inert ini Coupling 

1FSW0M = Residual Flexibility with Stiffness Coupling, without Residual Mass 
RFIWM = Residual Flexibility with Inertial Coupling, with Residual Mass 
RFIWOM = Residual Flexibility with Inertial Coupling, without Residual Mass 
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CHAPTER II; SHORT-CUT METHODOLOGIES 


ORIGINAL FAGE \S 
OF POOR QUALITY 


I . INTRODUCTION 


In Chapter I we discussed several modal coupling techniques. All these 
techniques necessitate the solution of the coupled booster/payload equations, 
i*e. they are **full-scale** methods. As we pointed out before, this solution 
is quite expensive^ especially if it has to bo repeated several times e.g. 
during a design effort. Although mass and stiffness changes during a design 
effort are often small, current practices used to design payload structures 
require a new ”full~8cale** solution every time such small changes in the 
payload are made. A similar situation exists in the case of payloads that an 
designed lor multiple flights with s»od<rrate coni igu rat ion changes. 

A need exists for the development of ''sfiort-cut** methods- Th»‘ tt‘rm 
'‘short-cut** method implies that the method should be able to evaluate small 
changes in the payload in a relatively short time. First, a short-cut method 
should avoid the direct solution of the coupled equations of the 
booster /payload system. Secondly, it should avoid as much as possible the 
interfacing between different organizations* This means that one should 
strife towards as much independence for the payload design organization as 
possible. 

The objective then of Chapter II is to present several of the most 
promising of these short-cut methods. Also an assessment of their strengths 
and weaknesses will be made. The first oc these methods will be discussed in 
the next section. 

2. THE PERTURBATION TECHNIQUE 


In this section we shall discuss a short-cut method which is based on a 
well known perturbation technique. We shall first disetiss the per I ui i»at J 
technique in general terras and then apply it to the particular probh^m of a 
booster /pay load system. 
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Let us consider a set of equations of motion of a certain structure, 

where [M^] and [Kq] are the mass and stiffness matrix of tne structure, 
respectively. The vectors |x} and|F| are the generalized discrete 
displacement and force vectors. The eigenvalue problem associated with Eq. 
( 100 ) can be WTi.tten as 

(-“O ['^o] ) {*01 ■ >''• 


( 100 ) 


The solution of this eigenvalue problem yields a modal matrix [*^03 and a 
diagonal eigenvalue matrix satisfying 


( 101 ) 


’ [‘J’ol [S]C‘**o3 " (J02) 

Sext, let us assume that the elements >f [w^] and [kJ undergo small 
changes, so that the new system of equations can be written as follow 


[m] { X ( + [k] }x} = {f} 


( 103 ) 


with 


[M] = [Mq] + c[MJ . [k] = [k^] + , [Kj] 


( 104 ) 



where ^ is a sAall parameter such that 


= [«i] - r.“] - [kq] • ' w * 


(105) 


Note that the matrix differences on the right-hand sides of Eqs. (105) are 
small) so that it is easy to determine a small H so that and [Kj] are 
of the same order of magnitude as [m], 


The objective of the perturbation technique is to obtain a solution for 
the new eigenvalue problem associated with Eq. (103) 

[m] + [k] j 

without actually solving Eq.(106), To this end, us write 


{x| - {xq} + e |x^} + {x^\ + ... 

HI = Hoi "" Hii 

2 

w * o)q + e + e u)^ + . . - 

This can be done because of the small changes in [H^J and [K^ J as 
expressed in Eq. (104)« Also, 


( 107 ) 

( 108 ) 
( 109 ) 






+ 



( 110 ) 


' .... 

OF POOR QUALi . 
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vhere I and I q I are the normal coordinate vectors of the unperturbed and 
perturbed system, respec tively, i.e. 


O ! i ) 

M ‘ [O ill ,1,2) 


with 


[♦I’ rMlE*] -to . [*]’'[ ■ t“'3 ,,,,, 

First, let us substitute Eq. (112) into Eq. (103) and preraultiply by [ 4 , ] * 


yielding, 


1,1 * - [* 3 ^ |f| 


(114' 


Substituting Eqs. (108), (109) and (UO) into Eo. (114) and equating 

coefficents of like powers of r , we can write, 

T 

Wo! * [“'nl Wol ■ [*o3^, <115) 

I’ll C"o3 I’ll • [*i3 wl - ^ [“oK“i,3 I’ol 

It is now possible to soWe Eqs* (115, 116, etc.) in sequence. The first 
Eq. (115) represents the unperturbed equation of motion, l.e., the modal form 
of Eq. (100). This solution is available or can be determined. 


(116) 


Once the vector |q [ is determined one can solve Eq. OlG) if and 

arc known. The determinat lort of these matrices is the suhjt'ct of next 
paragraph. 
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First, it is always possible to write *8 a linear combination of 

the eigenvectors , 

['t’J = [V> 

(117) 

because i-8 a complete set of vectors (i.e, they can be used as a basis 

for a vector apace). Mote that px] represents the coefficient matrix 

L\)J linear combination and must still he deterrained. to this 

end, let us introduce Eqs. (104), (108), (109) and (117) into Eqs.(113), and 

oniy keep terras in e and e , 

T 


{[♦o] * '[♦o] [“]) ([«o] + " [«i3) ([*o] ^ [i] 

k] ^ ' w ^ ' ['ll ) ( w ^ ' w w ) ,, , ~2 

■([«„]+ if-iJju 

. ^ A. • ^ A 1 ' _ - 1 ,1 . t • _ £ 1? 


(118) 
19) 


Equating coefficents of like powers in and b. we obtain from Eqs. 
(118-119), , T T 

W + [«] - - [to] ["il [Vl 

2 [.„] * [«]^[«^] + Kjcg 

Equations (120-121) can now be solve for • This erables us 

to solve Eq. (116) for | } , so that from Eq. (110) we obtain the first order 

approximation 


( 120 ) 

(121) 


|q| = 




( 122 ) 



Or 


QuALir/ 
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and from Eq. (112) we obtain] x|where we use 

[t] = [<t»j 

This perturbation technique can now be applied to a booster/pay loaa 
situation* The assumption is that only changes of order S are made in the 
payload, i.e. 

[Mp] - [%] + ' [Wpi] ■ - [‘S-0^ ^ " [SJ 


( 123 ) 


(124) 


where [Mp^] and mass and stiffness matrices of the 

unperturbed payload po. Let us write the mass and stiffness matrices in Eq. 

(28) for the perturbed payload P, 

T ■ T 


■"b * Wp, ^p"p'p 


with 


[Tp] = 




0 . Sp 
LO IJ 


1 ,T, 


IpMp^ 


_S_1!pVp 
0 : 


• [2p] 


= -[ 4 ] [ 4 ] 


(125) 


(126) 


Using Equation (124) we can write 


w • -(cO - ^ 


(127) 
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and, keeping terms in ( only, We obtain 





■ 

- (tC3 ■' 

' - ^ " cC] 

[«]) ( [<"] 

+ t: 

[<;]) 

1 (128) 

or 

- o " 

cO - e 

(to " 

[<;]- ccj " 

[■fi] CC3 

" (Cl) 

[Sp] 

l^po] * 






(179) 

with 


1 







' - to 






( 1 »() ) 



'' r<i3 



tC3 


(131) 

It J 

r. now possible to write [^|J n 

s 





= 1 
1 POJ 

^ [^PIJ 





(132) 


with 
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Finally we can write the mass and stiffness matrices of the coupled 
boos ter /pay load system as follows 

["cJ • ["a.] + ' ["c,] • ['‘c] ■ [\-oJ * ' 

where C indicates Coupled. Matrices and are given by expressions 

(125). Taking into accitunt Kqs. (124) and (1321 we ran write 


( 134 ) 






*^B *PO*S»O^PO ' 


> T 

'I K T 

PO PO 


1 


and 


(135) 



^TO^0^P1^^P0^1^P0''’^P1^0^P0 I ^^PO^l^^Pl'So^ ^PO 
_ 

- ^PO '^^O^Pl^^^'pi^PO^ 1 ^PO^l^PO 


(136) 



T?oWpi:^?oSiTpot!n_So^pol. 


1 ^PoSi’po 


(137) 


where [^po] and I^^Plj given by Eq- (133). Equation (134) is now 

equivalent to Eq. (104) and the perturbation technique can be applied. 


Note that theoretically one can also obtain the higher-order 
perturbations ^ 3 ^ etc. but for all practical purposes only • 

per lurbal i OUH are icu ludt»<l. The ilion is, how important t * , i \ 

etc. terms are. It is evident that Eq. (107) is only valid as long as the 
asymptotic expansion in f dors not break dowtii i.e, as lon^, as t |xj| 1 ^ol« 
‘ { ^7 f ’ l^|l» etc. 1’herr are indeed cases wlirre such an asvnipL®tlc exji.'mt; ( o.» 
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breaks down. It is then necessary to introduce oth^r perturbation techniques 
(e.g. Lighthill, multiple scales, etc). In this connnection it is important 
to roco^jnize the fact that small changes in the mass and stiffness of the 
payload produce small changes in the eigenvalues and eigenvectors, but not 
necessarily small changes in the response. These perturbation techniques show 
some promise and will be investigated further as part of Task II Methodology 
Devolopmont . 


3. THE BASE DRIVE TECHNIQUE 


In this section ve shall discuss the Base Drive Technique as developed by 
W. Holland, A. Devers and H.Harcrow. Let us first recall Eq. (87) in 
partitioned form, 


- L'Vpjr-St - 



(HK^ 


(I »y) 





I ' [■‘b"b‘'b*‘'p"p''p] ' ( - 


l‘'Vp] I -I I - ['X'p] I «t"'i 


i-'here we 


jolved for I in Eq. 


(140). 


The oayload designer is primarily interested in predicting the response of 


the pay'oad (see Chapter I, section 8.) which is given by 


4 with 



(141) 


— r , 


F 

1 s 


X.. 


N 


ny 



' y-. 


(142) 
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where 


-P 


BR| 


BF| 


jandjx^ j must be computed from, Eqs, (138-140) and from Eq. (86). 


The idea of a base drive short-cut method is to approximate 
(140) in such n way as to av<iid the solution of the complete se 
To evaluate a particular short-cut method, the approximation of 
compared to the exact value given by Eq. (140). 


/ ; in Eq. 

t ( 138-140). 
i 

I |must be 


A significant simplification of Eq. (140) occurs when the interface is 
statically determinate, when 

[^^b^b] ■ ['pVp] ' H 

This eliminates the dependence of 


(143) 


^rl 


and 


I I “^‘^nraes, 


1 


“BRl 

Xi f 


r T T T ^ 

IWb+^p"p’'pJ 


(-[^pVpjiift -[-pVpji^^i 
- b] I^n'I ) 


A first possibility is to assume that the presence of the payload has no 
effect on the response of the booster, i.e.jj^^^} = |o}*This approach is called 
the Direct Base Drive Technique. Indeed becomes 


[iH=p] 


[^Vp] 


Ixft 


(145) 


which means that the payload is "directly" driven at its base (i.e. its 
interface with the booster B) by the force on the right hand side of Eq. (145). 
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Theoretically, the D^^ect Base Drive Technique aasuaes there is no coupling at 
all between the booster and the payload* Practically, it means thatjx^^i << 
or that the feedback of the payload is negligible* The conditions under which 
such an approximation is valid is scill an unanswered question. The 
development of a criterion of validity of the use cf the Direct Base Drive 
Technique should be part of Task II: Methodology Development. This topic 
together with some other considerations will be discussed in Chapter III. 

4. THE IMPEDANCE TECHNIQUE [7l] 

In this section we shall discuss yet another approach to the solution of 
the equations of motion of the booster/payload system* The Impedance 
Technique as developed by K. Payne is basically a full-scale method in the 
sense that it does not make any assumptions concerning the size of the payload 
nor the extent of the changes made in the payload* However, the method does 
avoid a full-scale solution of the coupled booster/payload equations of motion 
and is particularly suited to deal with small changes in the payload. The 
Impedance Technique is essentially a Base Drive method (see section 3*). It 
differs from the approach in section 3* in the manner in which the interface 
accelerations I x^^jare computed* Indeed, the interface accelerations will be 
computed in the frequency domain instead of the discrete time domain thereby 
essentially converting a set of differential equations into a set of algebraic 
equations* 


Let us now derive the necessary equations* First, recall Eq* (3). 


T. 



«nd vriCe the top end bottom pertitions eeperetely, 

["p] IM ^ [s]wt ■ I- ■>?! 


(146) 


(147) 




(148) 
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where we invoked Eq. (2) and assumed the external payload forces to be absent 
(l^pl “ |of) • we int'^oduce an equation similar to Eq. (63) | 

* l*bl -"iV (14 

where the F vector represents the booster i.^sponse due to the externally 
applied force vectorj and therefore satisfies, 


° tfsl 


(150) 


and the R vector represents the response of the booster due to the feedback of 
the payload and satisfies 


[«b]|-b( '‘*[s])«b(' -{-o?! 

Applying transformation (33) to Eq. (150) yields, 



We now consider Eq. (149) and recall Eq. (4), 


(151) 


(152) 



( 153 ) 
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Keeping che bottom partition in Eq. (153) yields 



1*1 



B /• ^ 

Xj (t) 


R 


( 154 ) 


vhere ve now included the dependence of the vectors on the discrete time 
t. Taking the Laplace Transform (with zero initial conditions) of both sides 
of Eq. (154) aid introducing the transformation # “ jf^‘(withs* Laplace 
variable; '* the ith input frequency; and write 




\ B 




I B 
1^1 


(jn^) 


i = 1,2, 


(•55) 


(* number of input frequencie 


which represents Eq. (154) in the frequency domain. Taking the second time 
derivative of Eq. (155) yields 






Xt (jfi.) 


1 - 1 , 2 ,, 


(156) 


The basic idea of the Impedance Technique is to calculate the interface 

i 11Q.)lin the frenuenev domain and then transform them back to 


accelerations I (jQ^) jin the frequency domain and then transform 
the discrete time domain. The two terms on the righthand side of Eq. (136) 
will be replaced by algebraic matrix expressions so that the calculation 
of^- ® does not involve the solution c£ a set of differential equations. 

Let ua start with the first term on the right-hand sid^ of Eq. (136). To this 
end let us convert Eq. (l.*>2) to the frequency domain^ 


}’b «“,>}' 




|"b 


t - 1.2., 


(157) 
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or, taking' the second time derivative. 


1 — 
CM 

c: 

n ^ 1 ,1- 

-LrtJ 

1 D 

‘bJ fB^^Vj- 1 ^ 1 , 2 , 


Let us write Eq. ( 33 ) as follows, 



with 


r 1 


1 

__s_ 

♦b 


,B 

L J 


L J 


(158) 


(159) 


(160) 


Then the bottom partition of Eq. (159? reads, 


1*1 


I’B 


(161) 


I 



r 

M 

lU «( 


B I I B 1 I ’ 


1 - 1,2, (162) 


01 




[^A (JO^)j 


|f (jO ) j 

< ® M , 1 - 1,2,, 


(163) 
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f/1 



u 

♦l 

-^2 2 


\h\ 

L J 

■ “b 


L J 






(164) 


Equation (163) theni yields the first term on the right-hand side of £q* 
(156)* The TStrix k in Kq. (164) is the transfer admittance from the 

points of application of the external forcesj | to the interface 
accelerations. 


Similarily, Fq. (151) can be transformed into 


i •• ^ 

\\ Vjfi 



1 


B (j^^) 




(165) 


where this time 



(166) 


is the metric of coefficients for the point admittance for the booster at the 
interface. Equation (165) yields the second term on the right-hand side of 
Bq (.56). However, the reaction vector |rj (jJi^) | is not known a priori. To 
determine this vector let us consider Bq. (148) which represents the paylord 
equations of motion, and write it as 




"P 

+ 

p * p 

1 P 

-X 

» 3B 1 



1 


T 

«N 

1 P 




(167) 


where we used Bqs. (5) and (7). Introducing the modal transformation (40) and 
taking into account the properties (41) we obtain from Eq. (167) 
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wich 


(168) 


b^] ^ ['^ij 


(169) 


(170) 


and 



given by Eq. 


statically determinate* 


(30). 


Mote that 



is zero wran the interface is 


The top and bottom partitions of Eq. (168) can be %rritten as, 


• 


• [■;] l?l • [■:]!■: 


i - - •r‘’I 

I I M 


(171) 


(172) 


[4}[o]i 


We shall now asstone that the interface is statically determinate ( 
and calculate an expression | from Eqs. (171--172). First, we transform 

Eq< (171) to the frequency domain, 

F 


( 

1% 




o2 -2 
- Wp 


[«!] f 


x‘, OOl) j . 1 . 1,2,... 


(173) 


.nd from Bq. (172) we obtain, 
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[>^]’ [^']|% <i»i>j^['^j{«' ov| --|’^«v{ 


1 • 1 ? 


Substituting Eq* (173) into £q. (174) yields 






^2 -2 
“P 


51’Wf: ■».' 


+ |x^ av I = - ^^v|. i- 


from which wa obtain the following expression for /aq \| 

I I i'l 


{rj (jv| = [c <jv] j*i 


i ' 1,2, 


where 




I' o' 


o2 -2 

■ “P 


[<Y [<] - B] 


is the impedance matrix of the payload at the payload/booster interface. 
Finallyi we substitute Eq. (176) into Eq. (165), 



«”i> I 



Combining Eqs. (156), (163) and 

(H - [‘ t 


(J«i)j jxj (jn^) j ^ J 2, 

(178) yields 

1*1 0‘'l>{ * [* «"l>] {^B 


1 


1 , 2 ,, 


(174) 


(175) 

1 , 2 ,... 


(176) 


(177) 


(178) 


(179) 
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whert ve tlto ut«d Eq. (1) 

The coefficient oetrix of|*I in Eq. (179) represents the coupled 

iapedance of the booster /pay loed systea, and the ri^t hand side represents a 
pseudo generalized force. The interface acceleration can now be computed from 
Eq. (17*') with relative ease. 


If 

force I 


we now consider a new payload on the same booster and with the same 
Fg|, the ri^t hand side of Eq. (179) does not change so that, 


"B 


(jfi. 


:>l“- ([■]-[■'•• 



B (jJ2,) 


1 ( 1 ) 




( 1 ) 


provided the interface does not change. 



( 1 ) 




( 180 ) 


The interface accelerations in the discrete time domain can now be 
derived from Eq. (179) or Bq. (ISO), using the inverse Laplace Transform. The 
payload response then follows from Eq* (171). 


The approximation involved in the Impedance Technique is imbedded in the 
transformation to and from the frequency domain* If these transformations 
were exact, the method of determlniogj x^t^ould be exact. Therefore, one of 
the objectives of Task II* Methodology Development, shoul. * a detailed 
investigation of these transformations* There are also pr. .is pertaining to 
the modal damping when working in the frequency domain. 

Although Eqs. (179) and (180) were derived for an undamped statically 
determinate sy..tem, it is clear that damping and statically indeterminate 
interfaces can be included. For an indeterminate system it becomes necessary 
to keep track not only of the interface accelerations but also of the 
velocities and displacements at the interface* The use of the Fast Fourier 
Transform in obtaining the spectral data to be used in Eqs. (179-180) also 
presents some problems* In general, however, enough correlation with the 
exact time domain solution is apparent to warrant further investigation into 
po&sible improvements. 
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5. THE GEHERALKED SHOCK SPSCTEUM TKCHMIQinS ^2, 7^ 

The epproech preeented in Chit section is e generalised version of the 
shock spectrum technique as developed by Bamford C^23. 


In order to explain the basic ideas underlying this technique let us 
recall equation (44) 




--F 


y. 


* ♦rJVp'n " 


yj 


U\ 

4 




N 

0 


0 

(181) 


and let us assunc we reta.fned M nodes for the booster (ie is an MxM matrix) 
and N modes for the payload 6.e|^l is an NxN matrix) 

L 


The basic idea of the shock spectrum technique is to determine load maxima 
without having to solve Eq. (181). To reach this goal» a new model both for 
the coupled system and the forcing function is introduced* 


First, the (8 * M) modally coupled Eqs. (181) are replaced by (N x M) sets 
of two simultaneous equations each of which represents the coupling of one 
payload sxxle with one booster mode, as follows, 
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where we sssumed thst the interfsce is ststlcelly deteminsf 
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Secondly, s bound “lyp on esch of the (N x M) model reeponses of the 
psylosd is estsblishsd. This is done by introducing s new model for tht 
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forcing function in Eq. (182). Tue rather complicated forcing function is 
replaced by a much sioipler function (e.g. an impulsive force) which produces 
the same maximum response peak as the original force. Thereforei an 
analytical solution for both the response and maximum response of Eq. (182) is 
possible (after some additional simplifications). Finally, a bound on 

the total payload response Is constructed by summing over all the individual 
modal bounds vover absolute values or in a root-sum-square sense that 

allows for phase weighting). Payload member loads are obtained by adding the 
contributio ^ of all payload modal loads. 

As stated above, the forcing function in Eq. (182) is replaced by a modal 
delta function of a certain magnitude • This magnit*'de is evaluated 
from an already existing transient analysis of the booster with or without a 
dum^ payload. 

The main objections that can be raised against the Shock Spectrum Approach 

are: 


1) No critical evaluation is available on the validity of replacing 
model (181) by model (182). What effect does this replacement have on the 
load bounds? This change of modv^.l could not only result in a too conservative 
design but also in an unconservative one. Model (lb2) not only ignores the 
coupling between the B-modes due to rigid body feedback of P , but more 
importantly it ignores the effects that the coupling of one B-mode with one 
P-mode has on all the other F -modes. 

2) The manner in which is calculated again ^ ^aves the question of 
whether or not the envelope values are conseirvative or not and by how much. 

3) The c:.e.;hod appears rather complicated and is not simple to use. This 
can lead to misinterpretation and confusion when the method is applied. More 
rigor in the mathematical formulation is desirable. 

For these and some other minor reasons we do not favor further 
investigation of the Shock Spectrum Technique unless a better version appears 
that answers our basic objections. 


71 



CHAPTER III! PROPOSED IMPROVEMENTS OF EXISTING SHORT CUT 


METHODS - A NEW APPROACH 


1 . INTRODUCTION 


In Chapter I we reviewed and assessed four prominent "full-scale** 
methodologies* This allowed us to introduce the necessary background material 
in terms of a unified nomenclature* After careful evaluation we can state 
that all four methods have their merits* However t the Residual Mass and 
Stiffness Mechod appears to be the most efficient general purpose method, as 
discussed in section 9 of Chapter I. It is the full-scale method which best 
describes the booster structure in terms of a minimum number of modes, given a 
certain cut-off frequency for the externally applied force |^b| • fact 

that no payload information is required tc obtain the booster siodel is a very 
convenient feature in connection with the present study* Therefore, the same 
booster model can be used as long as che booster does not change* 


The Residual Mass and Stiffness Method can be used for comparison purposes 
when future short-cut methods are evaluated* Moreover, most short-cut methods 
require a full-scale ’•start-solution*’ before they can be applied, which 
demonstrates the need for an efficient full-scale method* Also, some of the 
short-cut methods (e.g* Base Drive Method) are based on their full-scale 
parent method* 


In Chapter II four general short-cut methods have been discussed and 
evaluated* Althougli each of these methods has its own merits, it is believed 
chat none of them is acceptable in their present stage of development to 
function as a standard short-cut method for general use* It is the purpose of 
this Chapter III to propose several possible improvements of these 
techniques* In additica, we also wish to present a new approach* Although 
still being developed, this new approach shows great promise* 


2* THE BASE MOTION TECHNIQUE 

The Coupled Base Motion Technique as explained in Chapter 1-7, leads to 
the fundamental set cf equations (138-140), which we repeat here for clarity 
of presentation, 
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Bv»] R'l ^ -- Rv„]r»fi 

*[4vp]i^i -- M(i»ri ^!?i) 

j»f 1 ■ ['Ivb^Vp] ' [^Jvp] j-rj - [■'pVp] hfl 
-^ipjl^'i -[^Kvwp] K’ 




As mentioned in Chapter 1^8 i the payload designer is primarily interested 
in the response of the payload i.e.^ 
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andjXp I must be computed from Eqs. (183-185) and |x^ jfrom £q. 

(86). The base drive method focuses on Eq* (184) which yields ! 

providedjx^ jon the ri^t hand side of Eq« (184) is known. The idea is to 

produce an expression forj^^^j without actually solving the coupled set of 

equations (183-185). 

I I «»pp I 

x^|Bnd|x° j in Eq. (185). These 
coefficients represent the ratio of the payload mass and the total vehicle 
mass. In many STS applications this ratio will be rather bjiall ( <10%). 
Therefore^ a first possibility is to ignore these terms in Eq. (185). 
Secondly, in many applications we can assume a statically determinate 
interface, ^ ®o that Eq. (185) becomes 
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(188) 


Ordinarily I the coefficient matrix of in Eq. (188) is not small and 

cannot be ignored. A first possibility is to assume thatj ^^jis small and can 
be ignored. This means that the feedback of the payload is not important. 

This can be a realistic assumption because the payload is usually small 
compared to the booster. In this case we can completely ignore | jin Eq. 
(184) and wite 




i"p I 
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(189) 


Equation (189) is now effectively decoupled from Eqs. (183) and (183). 
Physically, ignoring the feedback of the payload means that payload and 
booster are not modally coupled. The question is, when such an approach is 
valid. This is one of the points to be investigated in the future. 


A second possibility is to scale the vector l^'^jin Eq. (188). Indeed, let 
us assume a full-scale solution is available for some payload Now, some 
relatively small jl changes are made in the payload P]^ to generate payload P. The 
assumption now, is thatl^^jis not much different from i.e. 


I^BRi . \^BK\ 


(190) 


Equation (188) for payload fl can be written as 
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’•'JVb] - - kvB-^^n^pAi] {*r}^ 


(192) 


Taking into account Eq. (190) it follows from Eqs. (192) and (188) that 


{’'i'} - [’’bVb'^^JVp] ^ (|«) 

which yields a scaled value for|x^^lto be used in Eq. (184). Again, one 
should investigate when such an approach is valid. 


One way to improve the Direct Base Drive approach is to use a perturbation 
technique such as the one discussed in Chapter II-2. This could at the same 
time reveal when a Direct Base Drive is valid or not* Indeed, an asymptotic 

. (^p1 2 

expansion of . x^> could reveal the magnitudes of the terms in c , e etc. It 
should then be possible to decide if and when the aero-order term is 
sufficient or not to represent the response of the payload. 


Another possible route of investigation is given by the modal form of Eqs. 
(87), namely Eqs. (94), 
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Once {*^1 .iis knom in Eqs. (194) and (195) we have a decoupled ayatem of 
equations which can be easily solved. The idea is to start an iterative 
process with^q^ 1'(0> (i.e. Direct Base Drive) and solve aud^qj^| in 

Eqs. (194) and (195) which would be easy to do as nKsntioned above. These 

t **BR*f 

j^from Eq. (196). This 

value is then used in Eqs. (194-195) and the process is repeated until a 
satisfactory solution is obtained. Hopefully^ a fast convergence scheme can 
be developed. In the meantime, further research showed that no iteration process 
is necessary and a direct solution is possible. We shall discuss this approach 
in a future report. 


In conclusion, it can be said that the cet of Eqs. (183-185) or any 
equivalent set offer several possibilities for the development of an adequate 
short-cut method. The main problem is the development of criteria of validity 
for the several assumptions that must be made. 


3* THE I MPE DANCE TECHNIQUE 


As stated in Chapter II-4, the Impedance Technique is basically a Base 
Motion Technique. A set of suitable interface accelerations are given by Eq. 
(179), 



!■] - [■ (js^i)j |^c(j«^ ]Y [*«v ]i Fb0!1i)|^ 


(197) 


or, when the booster, booster forces and interface do not change we also have 
Eq. (180), 




( 2 ) 


[i]- [b(jv] j ' 


(198) 

O) 


The lapedance Technique Is an exact method, in the sense that the 
intsrface accelerations given by Eqs. (197>198) are exact, fhe approximation 
lies in the transformations from the discrete time domain to the frequency 
domain and back again. Indeed, problems were encountered with regard to the 
modeling of damping and the application of the Fast Fourier Transform. 
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The outstanding feature of the Impedance Technique is that no assumption 
is made as to the size of the payload. Also, the relative ease with which the 
interface accelerations are determined is a very attractive property of this 
technique. However, further research is needeJ. in the area of modeling the 
damping and in the area of conversion from the time domain to the frequency 
domain and vice versa. Also, we would like to investigate the possible 
application of a perturbation approach in case the payload and/or Che changes 
in the payload are small. 

4. A NEW TECHNIQUE 

In the course of our Investigation and evaluation of several short-cut 
methods it was noted that many methods involve assumptions and approximations 
leading to either doubtful or cumbersome resalts. In addition, it is often 
very difficult to assess the effects of those assumptions on the response and 
the loads of the booster/payload system. 

The basic problem is to somehow deal with the coupling effects between 
booster B and payload P without solving an eigenvalue problem pert 'lining to 
the coupled boostev/payloau system. This is a difficult problem Indeed* Each 
of the short-cut methods discussed in Chapter II addresses this problem in a 
different way. However, the proposed solutions invariably lead to cumbersome 
mathematics and program coding. This observation led us to the development of 
a more direct approach which we think shows great promise. This new approach 
is easy to understand and easy to implement. It is oased on the work of C. W. 
White and B. D. Maytum D4I . Although the theory is still being developed we 
already can present the basic philosophy. 

Let us recall Eq. (28) 

(199) 
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which represents the set of equations or motion of the coupled 
boor, ter /pay load system. It is now assumed that a cut-off frequency is defined 
based on a Fourier series expansion of . FurtLermore, we clso assume 

that e.g. the Residual Mass and Stiffness Method was used to construct the 
following set of modally coupled equations^ 


^ 



where 



( 200 ) 




( 201 ) 


and the cut-off frequency was used to determine the size of landf^l 
according to e.g. Chapter 1-5. In other words, the size of Eg. (200) ^3 
already much less than the size of Eg. (199). Due to e.g. the Residual Mass 
and Stiffness Method, the reduced Eq. (200) still represents an acceptable 
model for the coupled booster/payload system. 


The first step of the present approach is to solve the eigenvalue problem 
associat»>d with Eq. (200), namely 



T T ^ T T P 

lH-^BTpMpTp*3 ‘ 
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0 Pi 



( 202 ) 


yielding a set of modes [f} and a aet of frequenciea^^| aatisfying 
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(204) 


(205) 


where |u| ajLc* the new normal coordinates* bub^tituting transfoonation (203) 
into Eq* (200) and premultiplying ^ and using Eqs* (204-205) we obtain 

the uncoupled set of equations ^ 

* [“1 1 “( ■ [’] " 

The modal matrix J and thu frequency matrix represent the medal 
information of the coupled boos ter /pay load syaten* The idea now is to change 
the payload and calculate the changes in [ 4 /Jand In other words^ we use 

the full-scale solution of Eq* (200) as a "star t» solution"* This approach is 
taken in most short-cut methods and as such does not detract from the present 
approach. For example, this full-scale s* lution could be determined at the 
beginning of a design effort and would stay the same ^or a^^l subsequent design 
cycles of a particular payload* 

Let us now consider a new payload PI, with mass matrixfKp^^jand st*f^Xk .*?6 
matrix [ . This payload PI could be totally new or just a modif rf 

the nominal payload P, as long as we h ive the same degrees of freedom . r*. bclh 
payloads 7 and FI* For this new payload PI, we re^^ace Eq* (199) by 
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is the new displacement ' -ctor. 

Let us now writs Eq. (207) as follows 
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^ f^P**P^pj 

['‘bb] " [^nSi^pij - [^Vp] 
[Sp] ■ ['nSi'pi] - [‘^'V'p] 



(Mpj - «p) 




(213) 


(214) 


(215) 


Note that in case the interface is statically determinate 


if in addition the 


geometry of the payload is not changed the 
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(216) 


(217) 


(218) 


(219) 


( 220 ) 


Also, note that if no changes are made in the mass the right hand side, of 
Eqs. (213) and (216-219) become CoJ and dimilarly, if no stiffness changes are 
made we hrve from Eq. (208) thatE^pl* t^p^Jand consequently Eqs. (216-219) are 
valid while Eq (220) becomes [kpp]* j^Oj although the interface can still be 
statically indeterminate. 
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Next, let us define the following transforaation 



X 

After substituting Kq (221) into Eq. (210) and premultiplying by fAj we 
obtain the set of equations that now replaces the set £q. (200) 
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( 222 ) 


The next step is to define the transformation 




(223) 
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Tht. trantfomation (223) ia now aubatituted into Eq. (222) after which we 


prenultlply 


r 

by I 


and invoke properties (204-205), yielding 
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This Eq (224) replaces Eq (206) « For convenience, let us denote 
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(225) 


(226) 


Matrices [mJ and [kj represent wcie perturbations in the mass and stiffness 
matrices and fli^of system (206). At this point, several observations can 
be made. First, it should be noted that it is very possible that certain 
changes in the payload will only affect a limited number of modes and 
frequencies. This means that several columns in[>kj and corresponding 
elements in^i^will not change after the changes in the payload are made. 

This reduces the size of Eq. (224). Secondly, in solving the eigenvalue 
problem associated with Eq. (224) it is possible to use a Ray leigh-Ritz 
approach with as the estimated start modes. The smaller the changes in 
the payload the better estimate will be and the less iteration's will be 
necessary to produce the new modes and frequencies of the pertucbed booster 
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B/Paylo«d PI ay«tem« An «ven bett«r st«rting umt of modea could be the 
solution to the perturbed eigenvalue problem with ell off*diegonel terms equal 
to sero (this is equivalent to Che first term in a Taylor series expansion of 
the perturbed system modes and frequenc^^ Thirdly, ve wish to investigate the 
possibility of truncating modes In [^3 ecccrding to the initially defined 
cut-off fvequencv, If this possible Eq. (224) could be reduced in sire by 
approximately 50% compar^id to the already reduced system £q. (200). This 
reduction would be ia acidition to the one due to unaffected modes as mentioned 
above* However, this question must still be carefully investigated* Finally, 
it is also possible that the modes are group«id in subsets which show very 
Tittle or no coupling between each other. This means that the eigenvalue 
problem associated with Eq. (224) can be replaced by two or more smaller 
eigenvalue problems, vhich of course reduces the computation time. 

There are additional advant 'ges to this method: simplicity of use; 

accuracy of results (e.g. this method could even be used as a full-scale 
method); possibility of using engineering judgement and experience; the 
possibility to identify changes required to meet certain frequency 
requirements; the pcssibility to change branch frequencies to decouple the 
load problem leading to smaller eigenvalue problems, the potential for 
significant computational time savings. 

Note that it is also possible to solve eigenvalue problem Eq. (202) 
for the perturbed system, using the coupled modes of the unperturbed system 
as a first guess in a Rayleigh-Ritz type solution. 

Finally, we will investigate the possible combination of this technique 
with the Base Motion Approach. 

5^ CONCLUSION 


The purpose of this report is to define existing rvethodologies, evaluate 
their effectiveness In analysing dynamically coupled structural systems, and 
to define an approach where a "short-cut** methodology may be derived in Study 
Task II. This goal was defined within our proposal: 
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"Our approach to validation of existing "short-cut" methodologies will 
bo to disect the accepted state of the art mathematical description of 
a coupled payload/booster system to Identify tl ^ various interactive 
forces that arise as a result of that couplir^ The merits of each 
methodology will be judged upon the degree to vmich they represent these 
interactions and in relation to the costs (in elapsed time and computer 
dollars) required by the method". 

The previous chapters have definitized the problem by reviewing full 
scale methodologies (Chapter I), assessment oi strengths and weaknesses of 
defined "short-cut" methods (Uiapter II) and our proposed approaches to define 
an accurate and useable "short-cut" tae^hodology . 

The following critical coianents reflect detailed evaluations of the 
reviewed techniques. It should be noted only sensory mathematical studies 
were conipleted on these methods. During evaluation of derived techniques 
Study Task III, these methods will be reassessed with respect to the 
derived techni^^ues. 

PerturbatiOii Technigi e 

The driving principal of th i technique is that small changes in the 
mass and stiffness results in small changes in the modal characteristics. 
However, ^hls is not a sufficient mathematical premise to assure small 
changes in th^^ response. This method ;>hows some promise and will be 
Included within the rae^hodology development of Study Task II. 

Base Dri v e Techniqv.o 

Analytically the Base Drive Technique assumes no coupling hetween the 
payload and launch vehicle. The principal rhort com'^ng of this method is 
a lack of adequate definition that this "structural feedback" is negligible. 
The development of this criteria Is planned in Study Task II. 

Impedance Technique 

This technique is similar to a base drive method with the imknown 
criteria involving frequemv transformations. This technique will ue checkcnl 
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out through application to a perturbation technique in Study Task III, but 
will not be seriously considered during Study Task II. 

Generalized Shock Spectrum Technique 

The assumption that the total response for a single coupled mode can 
be related by coupling a single payload mode with an associate booster mode 
is questionable. Although the applied forces are replaced by modal delta 
functions. It is not clear to the au::hor that this would assure couscrviil i vc- 
ness. In fact, our opinion is that these model changes could result in 
either i too conservative or unconservative design. The application of a 
generalized force even though based on historical data is difr^-ult to assure 
proper conservativeness. Utilizing these forces may eliminate frequency 
dependence of the force but what direct influence this has on the 
dynamic response is unanswered. Additionally, due to the numerous mathe- 
matical assumptions, more rigor in the mathematical formulation is required. 
The generalized shock spectrum technique could be useful in sitiiations where 
weight considerations are not crucial (e.g., static buildings in an earth- 
quake analysis) or dynamic situations, ac a first approximation in conjunc- 
tion with a more sophisticated method. F : technique will be reassessed 
during Study Task III, but will not be 'leriously considered during Study 
Task II. 

Planned Activity 

The detailed study plan has been revised and is included as Figure 3. 
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FIGURE 3. STRUG rURAL DYNAMICS PAYLOAD LOADS ESTIMATE STUDY PLAN 
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